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ABSTRACT 


There exists a demand and supply of goods at various 
stations at different times. There are constraints like capacity 
of the tracks, number of empty wagons available. Along with the 
cost of transportation there exists a cost associated with delays. 
Railways have to route and schedule not only goods loaded wagons 
but also empty wagons, so as to minimize costs and delays. 

A mixed integer linear programming formulation has been 
given. Literature survey is mainly related to minimum cost flow 
problems in network and vehicle routing problem. We have 
developed a heuristic for a single good wagon scheduling problem. 

The methodology consists of dividing time horizon in the 
problem into a number of periods. Now in each time period 
scheduling of goods wagons and empty wagons is done. In case of 
goods wagon scheduling we try to solve a min-cost flow problem for 
a cumulative demand and supply in that period. Then routes are 
found fitting to the solution of min-cost flow, by simulation. 
From the scheduling of the goods wagons we come to know about 
requirements or availability of empty wagons. Allocation of empty 
wagons, is done by solving transportation problem with a criteria 
of minimum shortest distance. Simulation used in the routing and 
scheduling of empty wagons, is based on choosing next available 
shortest route. Procedure is repeated to cover all the periods in 
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the time horizon. 
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CHAPTER 1 
INTRODUCTION 


1*1 INDIAN RAILWAYS 

The first railway started in India in 185>. It ran 
over a route of 34 km long, from Bombay V.T. to Thane. Since then 
Indian Railway's network has grown enormously. Today. railways 
run over a route 70,000 km. long. It is the largest railway 
network in Asia and also fourth largest in world. Operations of 
railways have grown tremendously. It has emerged as a largest 
organization in service sector in India. In the world Indian 
Railways is second largest railway under a single management. To 
carry out the operations of railways, it is required to employ ZO 
lac employees, which is more than any other organization in India. 
Even though Railways, have enormously grown, both passenger and 
goods traffic has grown at a much higher rates straining the 
resources of railways. This has imposed many constraints. 

The Indian Railways haul 40 crore tonnes of goods in a 
year. Mainly it is the freight movement, that contributes to the 
profit of Railways. In fact Indian Railway bears a loss of more 
than Rs . 1300 crores on passenger transportation. Hence it would 
not be wrong to say, that goods transportation 


is of primary 
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importance to the railways because of the revenue it earns and 
profits it contributes. So also the scheduling of goods trains is 
the area of interest of our thesis work. 

1*2 CHARACTERISTICS OF THE OOODS SCHEDULING PROBLEM 

The passenger trains run according to the fixed time 
table. Even the new trains which are introduced are scheduled at 
time mostly according to the convenience of the passenger. Very 
few trains are introduced every year. As such, doing the schedul- 
ing technically, has very little scope. Contrary to this, goods 
trains do not have fixed time table, not even fixed routes, halts. 
The railways have to do entire scheduling (both day to day and 
long run) of goods trains. 

The Indian Railways is owned by Government of India. 
Hence, many a times it has to look after not only the movement of 
the goods, but also the distribution of the goods. It includes 
goods like food grain, petroleum products, fertilizers, cement, 
steel so on. Hence the nature of problem is some what tougher 
than the case in which customer would like to send goods from one 
place to another. This is because it is required to consider 
empty wagon scheduling and also time delays and accompanying cost 


1 ncreases . 
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1*3 THE OUTLINE OF THE PROBLEM 

In a typical transportation problem from operations 
research, we have warehouses which can supply certain quantity of 
goods. Also there are markets which have certain amount of demand 
to be satisfied. The cost of transporting from each warehouse to 
every market (cost matrix) is known. Similar to transportation 
problem, we have in our railway network supply nodes (stations) 
and demand nodes. We know for each arc in the network the cost of 
transporting both loaded wagon trains and also empty wagon trains. 
It should be noted that in our case cost of transportation can not 
be simply specified as cost of transportation from a supply node 
to demand node. This is because we have to consider not only the 
cost of transporting loaded goods but also that of empty wagons. 

Amongst the many constraints we have the important one 
being the capacity of the rail tracks. expressed in terms of 
number of wagons / trains that can be moved over a track in 
particular direction in unit time. This makes our problem 
partially is similar to a min-cost flow problem, having many 
source and sink nodes, instead of just one origin and one destina- 
tion node. However, the complexity of the real life problem of 
railways is of orders magnitude more than that of a min-cost flow 
problem. The railways have to deal with many types of goods which 
makes problem somewhat similar to the multi commodity min-cost 
flow problem. Further various types of goods can be transported 
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only in particular type of wagons. Each possible good type and 
wagon type combination having its own cost of transportation. The 
most important difficulty lies further when the concept of time is 
introduced. Every demand has a time deadline, i.e. the latest 
time by which demand has to be satisfied. So also supplies are 
available only after certain time and that supplies can be picked 
up any time after they become available. This concept of time is 
very important which is normally absent in of the typical trans- 
portation problem, travelling salesman problem, min-cost flow 
problem , vehicle routing problem. In our thesis work, we have 
dealt with the difficulties faced by the introduction of time 
associated with demand and supply. 

In our context, scheduling involves deciding the number 
of both goods wagons and empty wagons or trains moving from a 
origin node to a destination node, including the path along which 
they move, also deciding the starting time from the origin node, 
and finding the travelling time along the path, which involves 
specifying the idle time if any at any of the nodes on the path, 
on its way to the destination. Giving the delays in reaching the 
destination if any in reaching the destination. 

Me use the entire model developed, to analyze how 
introduction of more resources like wagons can reduce time delays. 


cost etc. 
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1*4^ ORGANIZATION OF THE THESIS 

Rest of the thesis is organized as follows. In the 
chapter 2, we give mathematical formulation of the problem. Here 
we have mixed integer linear programming formulation. This 
formulation reveals the complexity of the problem. In the chapter 
5, we have reviewed the literature on related topics. In the 
chapter 4, we begin by giving the reasons for choosing the 
particular methodologies, from among the various methodologies. 
This IS followed by complete presentation of the methodology we 
have adopted. In the chapter we take up a illustrative 
numerical problem. Then we present the computer results obtained, 
by running the program. The chapter 6, shows how the model 
developed can be used for deciding the critical investments like 
that of wagons. Finally, we present the conclusions derived from 
our results. Also we have suggested directions for future work in 


the last section. 



CHAPTER 2 

FORMULATION OF THE PROBLEM 


2» 1 INTRODUCTION 

In the previous chapter, the outline of the Railway 
Scheduling Problem was presented. In this chapter we present a 
mixed integer linear programming formulation of the problem. 

There are demands and supplies for different types of 
goods at different stations on the railway network at various 
times ^ We have to schedule the goods wagons and empty wagons so 
as to satisfy the dentrainds and supplies. Decisions regarding the 
movement of wagons and the type of goods they carry are to be 
taken. Loading and unloading schedules of wagons over a period of 
time are also to be decided. 

2* 2 FORMULATION 

2* 2* 1 Not-atrlons 

Decision Variables 

WUL , , No. of Wagons of type w that get unloaded at i with 

1 kwt 

goods of type k on day t. 

WLikwt wagons of type w that get loaded at i with good k 

on day t . 

WLM. ^ No. of loaded wagons of type w that start moving from i 
1 jkwt 
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WLI ,, . 
1 kwt 

WEM . . 

1 jwt 

WEI . . 

iwt 

1 j kw 

CE 

1 ] W 


WAV 

w 

‘^ikt 

®ikt 

^kw 

CD., 

1 k 

I 

T 

W 


to j carrying goods of type k on day t. 

No. of loaded wagons of type w that remain idle at i 
with goods k on day t. 

No. of empty wagons of type w that start moving from i 
to j on day t. 

No. of empty wagons of type w that remain idle at i on 
day t . 

Shortfall at node i, for goods of type k at time t. 

Constants CDat-aD 

Cost of carrying from i to j the goods of type k in the 
wagons of type w. 

Cost of carrying from i to j empty wagons of type w. 

Time (days) required to move wagons from i to j . 

Maximum no. of wagons of any type that can move between 
1 and j on any day . 

Total no. of wagons of type w that are available at the 
beginning and also in the entire planning horizon. 

Demand at node i, for goods of type k, at time t. 

Supply at node i, for goods of type k, at time t. 
Capacity of wagons of type w to carry goods of type k. 
Cost of delay at node i. for goods of' type k. 

No. of stations (nodes) on the railway network. 

Planning horizon in days. 

Number of different types of wagons. 
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2 . 2.2 

Objective 

Minimize 


A Linear Programming Formulation 
Function 


E E E E CL. . .<2 WLM. ^ + r V 

1 j k w Mkw ^ Ijkwt ^ ^ 


T. (CE 


/■r" uiCM 


E E CO <i: del > 

i k ^ t 


Constraints 

(1) Empty wagons balance 

Sum of number of empty wagons unloaded, idle on day t-land 
arriving on day t, a node i , on a day is equal to the number of 
wagons loaded, idle, left on day t at node i. Here it is assumed 
that loading / unloading takes one day. e.g. Nagon unloaded on the 
day t-1 can not leave on the same day either as empty or as loaded 
wagon but can leave on next day or later. 


E WUL.. + WEI ,, + E WEM. , 

^ ikw<t-l) iw<t-l) " JlW<t-tji) 

= E WL.. . + WEI. . + E WEM. . 
r' 1 kwt 1 wt “ 1 jwt 

k j 


for al 1 i = 1 . . I 
w = 1 . . W 
t = 1 . . T 


S imi larly , 

<2> Loaded Wagons Balance 

Sum of number of wagons loaded, idle and arriving at a node 
i, on a day is equal to the number of wagons unloaded, idle and 
left node i . 


EWL.. + EWLI , ,, ,v+EE WLM . , , . . . , 

ikw(t-l> " ikw(t-l> — " }ikw(t-tjv) 


= E WUL.. 4 . + E WLI .. . + E E 4 . 

7 ^ ikwt " ikwt , ijkwt 


WLM 


j k 


for all 1 = 1 . , I 
w = 1 . .W 
t = 1 . .T 
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(5) Satisfying demand and supply 

The cumulative capacity of wagons unloaded for satisfying 

demand at any node i, should be greater than the demand existing 

for all the types of goods, at any time, 
ti ti 


E WUL 

t = 0 


1 kwt ' 


'kw 


+ DEL . , , 
1 kt 


> 


E D 

t = 0 


ikt 


for demand 
for all ti = 1 

V 1 . k. 


T 


Similarly, the cumulative capacity of the wagons loaded at a 

node should not exceed the supplies available at that node i, for 

all types of goods, at any time, 
ti ti 


E WL 
t = 0 


i kwt 


'kw 


E S . . for supply 

t = 0 for all ti = 1 . . T 

V 1 . k. 


(4) Upper bound on wagons in transit 

E WEM ... + r WEM . . . + E E WLM . , ^ + E E WLM . , 
^ 1 jwt j iwt “ ijkwt Y jikwt 
w w k w k w 


< U 


for all links (joining i&j) 
for al 1 t 


We assume that loaded wagon which has to be unloaded, 
can not remain idle at a node (station). It has to be unloaded 
and kept idle. In other words the loaded wagon can remain idle at 
a node only if it is to be transported to some other node sooner 
or later and not if it is to be unloaded. Hence no. of unloaded 
wagons has to be less than the no. of wagons which arrive at a 
node for same type of wagon carrying same type of goods 


WUL.. . 

1 kwt 


< E WLM . ^ . . .. . 

y j 1 kw( t-tji ) 


for all k , w & t 
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(6) Total no. of wagons in the network is fixed 

E E ^ ^ ~ ^^'^w ^ at t = 0 

I k i (those at starting) 

2.3 OBSERVATIONS FROM THE FORMULATION 

In the last section we have given a mixed integer linear 
programming formulation for the problem of goods and wagon 
schedul i ng . 

It may be seen that we have included cost for not 
meeting the demand in time. The costs CD . could be appropriately 

1 K 

chosen to reduce the delays. However, if we try to solve the 
problem according to above formulation, then we find that even for 
a small network, the number of variables and number of constraints 
is very large. For example, if we want to solve a problem which 
involves a railway network of say 10 nodes and 20 arcs, only 4 
types of goods to be transported, in ? types of wagons and our 
time horizon is restricted to 100 discrete units of time on which 
wagon trains can be scheduled. alone requires 16,000 

variables to be defined, apart from the other decision variables 
to be defined and constant data to be given. The total number of 
constraints is around 25,000. This becomes extremely difficulty 
to solve even if we approximate all the variables as continuous 
variables rather than integer variables. This is well beyond the 
limits of the widely available linear programming packages. 
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Solomon [293 has shown that a problem of finding a 
feasible solution to a problem of a similar nature (but with 
lesser complexity) is a NP complete problem. This issue is 
discussed in the next chapter on literature survey. We limit 
ourselves to developing a heuristic method (which is computation- 
ally viable) which will generate reasonably good schedules. Hence 
the main aim of our work, is to develop algorithm and write a 
computer program for the same. The chapter 4 is devoted to 
explaining the methodology developed. 



CHAPTER 3 

LITERATURE SURVEY 


1 INTRODUCTION 

We begin the literature survey by studing the mincost 
flow problem, where the concept of time does not come in. Later, 
we study the work of various shades of Vehicle Routing Problem. 


3.2 MINIMUM COST FLOW PROBLEM CmCFP> 

The initial description of the transportation problem is 
due to Hitchcock <1941). Since then, there have been a number of 
papers concerning the computational aspects of MCFP. Surveys of 
earlier work were done by Charnes and Cooper < 1 9^1 ) , Ford and 
Fulkerson C143, Oantzig C93, Busacker and Saaty C43, Hu [183, and 
Golden and Magnanti C163. 

One of the most important and active area of interest 
today IS the basic simplex solution techniques which form the main 
tool of the work here. The approach is to use a basis tree to 
implement a procedure analogous to the primal simplex algorithm 
for the general linear programming problem. Studies have indica- 
ted that this is perhaps the most efficient way to solve pure 
minimum cost flow problems. The early precursors to this work 
were Dantzig (19^4) and Charnes and Cooper (1951) who applied the 
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simplex technique to the transportation problem. In the seventies 
several authors have utilized the special representation of the 
basis trees to increase markedly the computational efficiency of 
the primal simplex techniques. Langley et al.C20] uses the triple 
label method in a primal simplex algorithm for the capacitated 
transportation problem. Glover, Karney and Klingman [153 apply 
the triple label technique to the pure network problem. They have 
also shown the superiority of the primal simplex procedure for 
network flow programming problems which was attested by others. 

Another approach, which has also been incorporated in 
this work, as an alternate solution strategy is the incremental 
flow procedure. An incremental approach is discussed in Ford 
and Fulkerson C14]. Assad A. A C13 gives comprehensive survey of 
the literature dealing with multicommodity flow problems. A 
decomposition and partitioning approach is used to solve the 
min-cost flow problem using arc-chain formulation of the problem. 
The arc-chain formulation of the linear minimum cost flow problem 
suffers from the defect that the chains in the network should be 
enumerated in advance. These may be enormous in number and thus 
preclude solving a large network problem. An early paper of Ford 
and Fulkerson Cl 53 suggested that the chains be generated as 
neeeded according to a column generation techniques. Tomlin C523 
adapts the procedure in the context of the Dantzig - Wolfe 
decompos i t i on . 
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VEHICLE ROUTING PROBLEM 

The vehicle routing problem (VRP) involves the design of 
a set of minimum cost routes, originating and terminating at a 
central depot, for a fleet of vehicles which services a set of 
customers with known demands. Each customer is serviced exactly 
once and furthermore, all the customers must be assigned to 
vehicles such that the vehicle such that the vehicle capacities 
are not exceeded. The VRP area has been intensively studied in 
the literature, seeMagnant i C£5]. Christofides [^] has given 
optimisation algorithms and also heuristics for the basic VRP. WE 
refer to Solomon and Desrosiers C51] have described various shades 
of VRP with Time Window Constraints. 

3*3.1 Vehicle Routing Problem with Time Windows 

Solomon C513 has given mathematical formulation in which 
VRPTW and many other problem problem variations are derived from 
pickup and delivery problem (see section 5.5.2). This formulation 
appears in the Appendix A. In the VRP with time windows (VRPTW), 
the above issues have to be dealt with under the added complexity 
of allowable delivery times, or time windows, stemming from the 
fact that some customers impose service deadlines and earliest 
service time constraints. In these problems, the spatial aspect 
of routing is blended with the temporal aspect of scheduling, 
which must be performed to ensure the satisfaction of the time 
window constraints. The service of a customer, can begin 



within the latest time when the customer will permit the start of 
service. The times at which services begin are decision vari- 
abl es . 

It IS worth distinguishing now between hard and soft 
time windows. In the hard time window case, if a vehicle arrives 
at a customer too early, it will wait. In addition, due dates 
cannot be violated. In contrast, in the soft time window case, 
the time window constraint can be violated at a cost, in a manner 
similar to dualizing constraints in Lagrangian relaxation context. 

The vehicle fleet size and mix, may assumed homoge- 
neous. However, many of the methods proposed can be extended to 
consider a heterogeneous fleet size. Furthermore, the number of 
vehicles used can be free, i.e., the fleet size is determined 
simultaneously with the best set of routes and schedules rather 
than being fixed a priori. 

Additional complexities encountered in the VRPTW are 
length of route constraints arising from depot time window. It 
may be noted that the vehicle departure times from the depot are 
decision variables. Furthermore, precedence relationships among 
certain customer are encountered in some problems, such as those 
involving both pickup and deliveries. 

Most of the effort has been directed at the operational 
problem of determining the best set of routes and schedules. In 
the presence of time window, the total routing and scheduling 
costs include not only the total travel distance and time costs 
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considered for routing problems, but also the cost of waiting 
time which is incurred when a vehicle arrives too early at a 
customer location or when the vehicle is loaded or unloaded. 
Several authors have also have analyzed the strategic question of 
the optimal fleet size. Keaton [21] dealt with designing optimal 
railroad operating plans. Various related decisions are modeled 
as mixed integer linear program. Solutions are obtained using 
Lagrangian Relaxation and Heuristic approaches. 

The VRPTW has emerged as an important area for progress 
in handling realistic complications and generalizations of the 
basic routing model (Schrage [27]). In the last few years we have 
witnessed the development of a fast growing body of research 
focused on the vehicle routing and scheduling problem structures 
with time windows. Given the intrinsic difficulty of this problem 
class, the later work on the generic VRPTW has typically focused 
on heuristic methods. The findings of Solomon [29] indicate that 
this problem class is fundamentally more difficult than the VRP. 
The VRPTW is NP-hard (by reduction from the VRP) and heuristics 
have so far offered the most promise for solving realistic size 
problems . 

Solomon [29, ?0] has designed and analyzed a variety of 
route construction heuristics for the VRPTW. The results of the 
extensive computational study reported in Solomon [29] indicate 
that a sequential time-space insertion algorithm proved to be 
very successful in a number of important VRPTW environments. Very 
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good initial solutions to a large number of 100 customer test 
problems were obtained. The excellent performance of such 
heuristics can be explained by realizing that while routing 
problems seem to be driven by the assignment of customers to 
vehicles as indicated by the success of the Fisher and Jaikumar 
[12], generalized assignment heuristic- the sequencing aspect of 
the problem seems to drive routing problems dominated by time 
windows. It is this aspect of the problem that insertion approa- 
ches capture so well. A very large computational requirement was 
also reported by Cook and Russell [8], in early work on the VRPTW 
for a k-optimal improvement heuristic. which was effective in 
solving an actual problem with a few time constrained customers. 
Analytical results concerning the behaviour of VRPTW approximation 
methods are derived through worst-case analysis of heuristics in 
Solomon [293. For a variety of heuristics, including improvement 
methods, it is shown that their worst-case behaviour on n-customer 
problems, is at least order of n for minimizing the number of 
vehicles used, the total distance travelled and the total schedule 
time. While heuristics have been found to be very effective and 
efficient in solving a wide range of practical size VRPTW. few 
optimal approaches have been reported, see Kolen et al. [22]. The 
former authors extend the shortest q-path relaxation algorithms of 
Christofides et al. [6] to the problem with time window. The 
largest problem solved to optimality involved 4 vehicles servicing 
14 customers with tight time windows. 
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Ferland and Fortin C113 have given a heuristic approach 
for the vehicle scheduling problem with sliding time windows The 
advantage has been taken of the fact that the starting time of 
each task must fall within a given time interval rather than being 
fixed. The approach consists of identifying pairs of tasks 
offering good opportunity costs for reducing the overall cost and 
finding the ways of modifying starting times to allow them to be 
1 inked. 

3.3.2 Pickup and Dolivary Problems CPDPTW3 

PDPTW is a generalization of the VRPTW which is concern- 
ed with the construction of optimal routes to satisfy transporta- 
tion requests, each requiring both pick and delivery under 
precedence, capacity and time window constraints. Note that the 
VRPTW IS the particular case of PDPTW where the destinations are 
all the common depot. 

Most of the literature on pick-up and delivery problems 
has appeared for the dial-a-ride problem <DARP> . It can be seen 
that the single vehicle DARPTW is a constrained version of the 
classical travelling salesman problem. Psaraftis C24] presents 
two dynamic programming algorithms using respectively backward and 
forward recursions. The objective function minimizes the total 
customer inconvenience and "maximum position shifts” define time 
window constraints. These algorithms require 0(n^j”) time (where 
n is the number of customers), a fact which limits the tractable 


problem size to no more than 8-10 customers. 
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Sexton and Bodin C283 give mixed integer nonlinear 
programming formulation for the single vehicle dial-a-ride problem 
with given maximum delivery times or specified minimum pickup 
times to minimize customer inconvenience. Benders' decomposition 
is used to partition the problem into a routing component (a 
linear program). For the scheduling problem, an optimal algorithm 
which exploits the network flow structure of its dual is develop- 
ed. The overall procedure is heuristic due to the difficulty of 
solving the routing problem exactly. Its computational tractabi- 
lity stems from the non- i terat i ve nature of the scheduling. 
Successful computational experience on moderately sized real data 
is reported. Problems with 7 to 20 customers were solved in an 
average of 18 seconds of UNIVAC 1100/81A CPU time. 

Jaw et al . C19] also use a parallel insertion algorithm 
which appears to be very effective and efficient in minimizing a 
weighted combination of customer dis-utility and system cost. The 
authors report computational experience with simulated and real 
data on the VAX ll/7?0. The former problems. involving 250 
customers and 10-14 vehicles, took around 20 seconds each. 

The Bodin and Sexton [53 procedure is a traditional 
"cluster first, route second" approach. For a fixed fleet size, 
it partitions the set of requests into vehicle clusters and solves 
the resulting single vehicle dial-a-ride problems using the 
heuristic based on Benders decomposition (Sexton and Bodin C283). 
Requests are then moved one at a time from one vehicle to another 
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while attempting to reduce total user inconvenience. 

3.3.3 The Multi -Period Vehicle Routing Problem CMPVRP3 

Here, the time windows are full days and a service 
activity must occur on a specified number of days of the planning 
VRPTW where each customer has to be visited a given number of 
times within its multiple time windows. Most of the existing 
literature is based on application that require the solution of 
large scale problems. 

Beltrami and Bodin C2] address this problem in the 
context of routing hoist compactor trucks. The problem, which 
involves customers requiring service either 5 or 6 times per week, 
IS solved heur i st 1 cal 1 y . One approach involves developing the 
routes first and then assigning to days of the week. In the 
second approach, routing is performed after the customers have 
been randomly pre assigned to days of the week. 

Rusell and Igo C263 extened several routing heuristics 
to assign accounts to days of the week such that the weekly travel 
distance is minimized. The refuse collection problem considered 
involves the routing 4 trucks through 490 location requiring 
service from 1 to 6 times per week. The original problem is 
transformed into a pure routing problem by creating multiple 
copies of each account that requires service several times a week, 
resulting in a 776 customer problem. The best approach found was 
to assign the locations to days of the week using a simple 
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heuristic, solve the resulting problem, where the spacing condi- 
tions are enforced and improve the solution through the use of a k 
interchange heuristic <M-Tour, Russell C2?]). The use of M-Tour 
proved once again to be computationally prohibitive. Christofides 
and Beasley C73 develop new heuristics for the MPVRP based on 
median problem and traveling salesman problem relaxations. The 
algorithms are characterized by the use of an interchange proce- 
dure which improves on initial choices of delivery days combina- 
tions for customers. 

Haghani A. E. C17] presents the formulation and solution 
of a combined tram routing and makeup and empty car distribution 
model. The model results in a large scale mixed integer programm- 
ing problem with nonlinear objective function. A heuristic 
decomposition technique was developed to solve the model. 



CHAPTER 4- 
METHODOLOGY 


1 INTRODUCTION 

The problem as formulated in chapter 2 is fairly complex 
and hence we begin solving it by relaxing a class of constraints 
at a time and solve the relaxed problem by known methods which are 
efficient. We present the solution method for a single commodity 
problem, but it can be extended to the multicommodity case. 

The detailed steps in the algorithm for solving Railway Wagon 
Scheduling Problem are as given below. 

4. 2 INPUT 

(i) Network is specified in the form of arcs. 

Data about the Arc 

(a> Nodes : nodel , node2 which are joined by the arc 
(b) Distance. 

Cost for carrying loaded wagons cost for carrying empty 
wagons is assumed proportional to distance and both 
costs are computed. 

Cost is expressed as cost/rake for travelling arc. 


Assumed to be same in both directions. 
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(ii> Demand / Supply of goods (expressed in wagons (trams)) 

Input : Node, demand / supply, time . 

(ill) Length of time horizon 

length of one time period within a time horizon. 

We partition the problem into a few distinct time periods 
say equal to *n* . This can be done by a suitable heuristic 
or by simple inspection of data. We have done it by 
inspection only. 

(iv) Availability of empty wagons at various nodes in time period 
one (or in any other period) 

Input : Node, availability of wagons, time. 

4*3 OENERATE THE INFORMATION FOR GIVEN NETWORK 

For each arc, there is similar arc in reverse direction. 
Hence, ours is undirect ional network for which we have following 
relations holding true. 

Number of arcs leaving a node i 

= Number of arcs arriving at a node i. 

= Degree of a node i / 2 

= difference between pointCiD and pointCi+1]. 

(1) Arc numbers of all the arcs leaving a particular node i are 
stored in the form of a forward star network representation. 

Arcs' leaving node i will be stored in position from pointCi] to 
pointCi+13-1 in the array ftraceCD. 
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<2) Arc numbers of all the arcs arriving a particular node i 
stored in the form of a reverse star network representation. 
Arcs entering node i will be stored in position from pointCi] 
pointCi+lD-1 in the array rtraceC]- 


4.4 SCHEDULING OF GOODS WAGONS 

Out of these *n' distinct period problem, we pick up 
first period problem and totally drop the time constraints 
supply and demand. The min-cost flow problem PI is formulated 
below. 

Minimize EEC. ..X. . 

^ ^ 1 j 1 ) 

(1) Flow conservation at all nodes i 

Outflow - Inflow = Supply (Demand) 

<2) Capacity restriction on flows over the arcs 

X . . < K . V arcs 

1 J 1 J 


(?) Non negativity of flow constraint 

X. . >0 V ( i . j ) 

1 ) 

Where 

X. . Flow from node i to node j 
1 J 

C Cost of transportation unit quantity from node 
) 1 

k. 

B. Supply <or Demand) at node i 

Capacity of arc from node i to node j 


to node j 


are 

to 


the 

on 

as 



2 ? 


The solution to the above problem gives the flows over 
the arcs in a period. Various fractions of the flow over the same 
arc can be a part of the different routes from various supply 
nodes to different demand nodes. Solve the min-cost flow problem 
as given below. 

4-* 41 Solv© t-h© min— cost- flow for t-h© curr©nt. period 

Find the cumulative demand / supply of goods within the 
current time period. This is done because demand / supply may 
occur at different times in same time period. The demands must be 
satisfied within the same time period only and also at time 
mentioned. However supplies not exhausted during past time 
periods are considered available in the current time period also. 

Cumulative demand = J] Demand for current period. 

Cumulative supply = 2 Supply for current period + 

residual supply of previous periods 
With this information about cumulative demand / supply at each 
node and the data about the network, solve the min-cost flow 
problem for the current time period. The capacities of the arcs 
used in solving min-cost flow problem are specified for entire 
period. The min-cost flow is formulated as a linear program 
and solved by a revised simplex method. 

The output of the min-cost flow problem is in the form 
of flows over various arcs. The flow may exist in both the 
directions on a arc (arcs have been separately identified for each 


direction). Since the individual arc flows have already been 
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decided the cost of transportation is same irrespective of the 
route chosen and the decisions about which supply node supplies 
towhich demand node, since it fits in solution given by the 
min-cost flow. 

At At 2 Solving t-he Transporiat-lon Problem for goods wagons 

Take a demand node for all those nodes having a demand 
and each time it has a demand. Similarly have a supply node for 
all those nodes having a supply and each time it has a supply. 
Thus more than one copies of same node will exist if and only if 
demand / supply exists more than once till end of current time 
peri od . 

Now we wish to know which node is to supply which demand 
node and also the likely quantity supplied. While doing so it is 
desirable to minimize the delays. We solve the following problem 
PJ which is the transformed transportation problem. 

Formulation of transportation problem is given here. 

Problem P5 

Maximize EES.. X . 

^ ^ 1 j 1 j 

Subject to 

<1> Demand restriction for node i 

£;Xj.>A. Vi=l,2,...m 

<2) Supply restriction from node j 

J] Xj ^ < B. V j = 1 .2 n 

<5) Non negativity of flow constraint 

X.. > 0 V <1.J) 



Compute a shortest oListcmjce matrix 
between each pair of supply node 'i 
it into a shortest time matrix <t 


(d ) giving shortest 
I J 

' and demand node ’ j ’ . 
> . 


di stance 
Convert 


I J 

Compute a available time matrix giving difference of time 

at which supply becomes available and the time at which demand 

needs to be satisfied. 

Compute a slackness matrix. (S ) 

1 J 

Slackness = Available time - time for travelling shortest distance 


Use the criteria of maximimum of slackness to decide allocation 
between demand and supply nodes. This criteria is used so as to 
minimize time infeasibi 1 ities in current and future periods. A 
simple heuristic approach is used to solve this transportation 
probl em. 

Solution procedure to solve objective function by maximum 
slackness criteria is given below. 

Heuristic One 

Step 1 Find out slackness matrix (also known as transportation 
table) 

Step 2 Select a variable with maximum slackness as basic 

variable. It is assigned a value equal to minimum of 
supply or demand. Update the residual demand and supply 
by subtracting the assigned value. 

If there is tie in the selection of the variable with a 


Step 5 
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maximum slackness, th&n find the penalties for not 
making that variable as a basic variable. Break the tie 
by choosing the variable which has maximum penalty. 

Step 4 Repeat the procedure till all demand is satisfied. 

We try to stick to the above solution of transportation 
problem P5 and find the routes from supply node to demand node 
within framework of solution given by min-cost flow. 


If we get a infeasible solution to the min-cost flow 
then we need to solve a different transportation problem which 
will reasonabally take care of the costs of transportation. We 
solve a transportation problem P4. which minimizes the sum of the 
shortest distances. 

Formulation of transportation problem which minimizes the 
sum of shortest distances is given here. 

Problem P4 

Minimize E E D X ■ 

la L, j j , j 

Subject to 

<1> Demand restriction for node i 


E X. . > A. V i = 1 ,2, . . .m 

(2) Supply restriction from node j 

EX.<B. Vj=1.2....n 

<5> Non negativity of flow constraint 

X. >0 V ( i . J ) 

1 ) 

Where a shoTtest distanco matrix (D^^) giving shortest 


d i stance 
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between each pair of supply node ’i’ and demand node 'j'. 

Use the criteria of minimum of shortest distance to decide 
allocation between demand and supply nodes. A simple heuristic 
approach is used to solve this transportation problem. 

Solution procedure to solve objective function by minimum 
shortest distance criteria. 

Heuristic Two 

Step 1 Find out shortest distance matrix C transportation 
table > 

Step 2 Select a variable with minimum shortest distance as 
basic variable. It is assigned a value equal to minimum 
of supply or demand. Update the residual demand and 
supply by subtracting the assigned value. 

Step 5 If there is tie in the selection of the variable with a 
minimum shortest distance, then find the penalties for 
not making that variable as a basic variable. Break the 
tie by choosing the variable which has maximum penalty. 
Step 4 Repeat the procedure till all demand is satisfied. 

‘4.4*3 Simulating to find the routes for goods wagon scheduling 

We start from a supply node and travel along any of the 
arcs leaving the node if that arc has a flow indicated by min-cost 
solution. We recursively continue search till we reach a demand 
node which has some positive allocation for that supply node and 
demand pair. When a path is found successfully, we send a flow 
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along the entire path. Quantity of flow that can be sent is 
minimum of 1. minimum arc capacity along the path 2. minimum of 
supply and demand. In case we fail to reach such a demand node, 
then we backtrack to the previous node along the route so as to 
try with some other arc leaving the previous node. This back- 
tracking to the previous node along the route may have to be done 
recursively till we do not have any node left (other than supply 
node itself) in the route we are trying. In that case we give up 
the current supply node, which probably may have surplus supplies 
which may remain even after satisfying all demand. Alternatively 
its supplies are available too late and can not satisfy demand 
without delays hence considered later. 

The above procedure is repeated so as to find routes 
from other supply node to demand nodes. If the demand remains 
unsatisfied even after trying with all the supply nodes once, then 
we reject the tranportat i on solution which we developed earlier. 
It may be noted here that transportation solution developed 
earlier with criteria of maximizing the slackness was just to 
minimize any time i nf eas i bi 1 i ty that may occur during goods wagon 
scheduling or consequent empty wagon scheduling. Hence even if we 


do not stick 

to 

the allocation, it 

1 s 

always 

possible to 

find 

routes from 

the 

supply node to the 

demand 

nodes within 

the 

framework of 

the 

min-cost solution 

1 . 

e. without exceeding 

the 

flows indicated 

by min-cost solution 

which, in 

turn do not exceed 


the capacity of the arc specified for the period. 



If our problem is too tight problem then the infeasibi- 
lity IS found while solving the min-cost flow itself. In such 
cases we use minimizing the sum of the shortest distance between 
the supply and demand node pair as the criteria for deciding the 
allocations. Now the shortest available routes are found by 
simulation. Then, this approach becomes very much similar to 
empty wagon scheduling described in details later on. However in 
case of goods wagon, their scheduling is always done so as to 
reach the destination just in time (whenever possible). On the 
other hand in case of empty wagon, scheduling is done so as to 
start from the source node as soon as they become available. 

When sending the flow along a route found 

a. Check time feasibility if any, and find the amount of delays 

when it occurs. 

b. Keep the account of residual capacity of each arc on each day. 

Goods Flows are schedul&dt so as to reach dietnanxi node 
j'xist in time of need, rather than sending the flow at a time when 
supplies are auailaJble. 

This reduces the tightness in scheduling and time 
infeasibility that might be encountered during empty wagon 
scheduling. The exception to this scheduling rule. is the 
situation when supplies are not available within required time. 
Then the goods trains are scheduled to leave the supply node only 
when supplies become available, even though they will reach demand 
node with delay. 
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EMPTY WAGON SCHEDULING 


The demand and supply of goods (goods trains) acts as a 
source of supply and demand for empty wagons respectively. 
However the time of demand / supply need not be same as that in 
data. Delays might have occurred while satisfying demand of goods 
(which leads to time infeas ibi 1 i ty) . It is possible to have a 
situation in which demand gets partially satisfied in time and 
partially with delays. Also supply of goods may be picked up at a 
later time, rather than the time at which it becomes available. 
So also a smaller portion of the total supplies available at a 
node may be picked up at different times. 

Thus to get precise information for requirement and 
availability of empty wagons, it is necessary to keep account of 
(1) Number of empty wagons and (2) time at which they are required 
or are available as soon as loading or unloading of goods takes 
place during goods wagon scheduling. 

4.5.1 Solving the tranportation problem for Empty wagons 

Similar to goods wagon scheduling we shall form demand 
nodes and supply nodes for empty wagons each time there is 
requirement for or availability of empty wagons. Then we form a 
transportation problem to decide which supply node will supply 
which demand node. However unlike previous case the basis for 
allocations is simply shortest distance matrix which gives the 
shortest distance between every pair of supply and demand node. 
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The reason why we have used criteria of maximizing the 
slackness while goods wagon scheduling and that of minimizing the 
sum of shortest distance while finding allocations during empty 
wagon scheduling has to be noted. During goods wagon scheduling 
min-cost flow was solved first and the purpose of solving the 
transportation problem with maximum slackness was to minimize the 
time 1 nf eas i b 1 1 1 t 1 es that may occur during goods wagon and empty 
scheduling. Overall cost factor was already taken care of when we 
solved min-cost flow problem. The reason why shortest distance 
matrix is used as a basis for allocations during empty wagon 
scheduling is that we would like to minimize the costs of trans- 
portation for empty wagon scheduling also. We have also found 
that maximizing the slackness is a poor criteria during empty 
scheduling. Also it is not needed because nothing further depends 
on empty wagon scheduling, unlike the case of goods wagon schedul- 
ing from which empty wagon scheduling is derived. 

4.5*2 Deciding the routes for empty wagons by simulation 

Procedure for deciding the routes from supply node to 
demand node for empty wagons is somewhat simpler in case of empty 
wagon scheduling than goods wagon scheduling. Here we strictly 
stick to the allocations. Moreover unlike the min-cost flow of 
the goods problem, the capacities of arc are not expressed for 
whole period but for a each unit time within a period, which 
IS a fraction (to be specified) of the capacity for the period. 



We take a supply node. We try to send empty wagons 
demand node along a shortest path. In case, we find that there 
no capacity available on a particular arc in the path on 
particular day then 

(1> We try, keeping the empty wagons idle just before the arc 
which there is insufficient capacity available on that particul 
day till, sufficient capacity becomes available. 

<2) In case we find that by keeping the empty wagons idle, be 
that particular arc, it is not possible to reach the destinat 
in time, then we try to find a next shortest path excluding 
arc(s> whose capacity has saturated. In case we encounter simil 
problem in the next shortest path also then we repeat the proc 
dure. If the distance involved in travelling by next shorte 
available path is such that it will not be possible to rea 
destination then we accept that as time infeasibility, 
wagons are scheduled along a first shortest path only. We note 
delays that are bound to occur when we send the goods by 
shortest path only. 

The steps 4.4 and 4.5 are repeated for all the time periods in 
time horizon. The concise picture of the algorithm appears 


figure 4.1 




PIG. -1.1 CDIICIGE riCTURK OF ALGOHI'IHH 















36 


o OBSERVATIONS 

The choice of length of the time period seems to be 
important. If the period length is small then the number of 
periods in the horizon increases. and there is corresponding 
increase in the computational time. If the period length is taken 
to be very large then the wagons remain idle in a period. This is 
because wagons are scheduled only twice within a time period. 
Once goods carrying wagons are scheduled and once empty wagons are 
scheduled, so as to meet the requirements for the goods wagon. 
Hence there is poor utilization of the wagons when period length 
IS large. Also it may result in quite suboptimal and at times 
infeasible solutions. In which case we have try again with a 
smaller length of the period. 



CHAPTER 5 
RESULTS 


5.1 INTRODUCTION 

In the previous chapter, we have presented the solution 
methodology explaining the details of various steps involved in 
the algorithm coded. The program code has been attached in the 
appendix B. In this chapter we take up a numerical example. In 

the section 5.2, the data of the problem is given. Figure 5.1 
shows the railway network which is used for the problem. In the 
section 5.5, the output results have been tabulated. In the last 
section, we have discussion on the results we have obtained. 

5.2 NUMERICAL PROBLEM 

5.2.1 DaLa about arcs of the railway network 


Table 5. 1 Arcs of the network 


Arc 

No . 

Node! 

NodeZ 

D 1 stance 

< km> 

Capac i ty 

Problems 

1 2 

1 

Delhi (1) 

Jhans 1 

(2) 

400 

16 

26 


Delhi <1> 

Kanpur 

(5) 

400 

20 

25 


Kanpur <5> 

Jhans i 

(4> 

200 

5 

5 


Ahmad .(15) 

Jhans i 

(2) 

800 

5 

5 

5 

Kanpur (5) 

Lucknow 

(4) 

100 

25 

55 
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6 

Jhans 1 (2) 

Bhopal (19) 

500 

50 

55 

7 

Ahmad . < 1 5 ) 

Bhopal (19) 

550 

6 

6 

8 

Bhopal (19) 

Indore (20) 

200 

25 

25 

9 

Lucknow <4> 

Calcutta (5) 

1000 

15 

25 

10 

Bhopal <19) 

Bhusaval (15) 

400 

7 

8 

11 

Bhopal (19) 

Nagpur (16) 

550 

9 

7 

12 

Calcut. (5) 

Jamshed • ( 6 ) 

250 

15 

25 

15 

Bombay (12) 

Manmad (14) 

250 

25 

19 

14 

Manmad (14) 

Bhusav. (15) 

200 

16 

16 

15 

Bhusa. (15) 

Nagpur (16) 

400 

7 

7 

16 

Nagpur (16) 

Jamshed. (6) 

250 

25 

25 

17 

Bombay (12) 

Pune (11) 

200 

6 

6 

18 

Manmad (14) 

Daund (10) 

250 

7 

7 

19 

Nagpur (16) 

Balhar. (17) 

200 

19 

16 

20 

Oaund (10) 

Bangal . ( 9 ) 

900 

6 

6 

21 

Balhar. (17) 

Hydrab. (18) 

550 

9 

9 

22 

Jamshed . ( 6 ) 

Vishakh. (7) 

800 

6 

6 

25 

Banga 1.(9) 

Hydrab. (18) 

700 

6 

6 

24 

Bangal . ( 9 > 

Madras (8) 

550 

25 

25 

25 

Hydrab .(18) 

Vishakh. (7) 

600 

25 

21 

26 

Hydrab. (18) 

Madras (8) 

700 

5 

5 

27 

Madras (8) 

Vishakh. (7) 

850 

5 

5 

28 

Pune (11) 

Daund (10) 

100 

5 

5 

29 

Ahmad . (15) 

Bombay (12) 

500 

25 

25 



I, j : nodes 
r : arc number 


RAILWAY NETWORK USED IN THE PROBLEM . 
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Network shown contains 20 nodes <n) and 29 arcs (m> . 
There are also 29 more arcs in direction just opposite to the one 
shown in the network. These are identified by numbers from m+1 to 
2m. The cost of the wagon movement along a arc is taken propo>r- 
tionai to the distance. 

The cost of loaded wagon train movement = distance (km> / 10. 

The cost of empty wagon train movement = distance (km) /l^. 

Also to have the consistancy in the uni tsfol lowing conversions is 
used . 

1 unit time = 4 hrs = 200 km. 

1 day = 6 units of time. 


5* 2* 2 Dat.a about th© availability of supply 

Table 5.2 Supply at various nodes 


Supply node 

Supply, TimeCday) 

2 

Jhans i 

9,2; 10,3; 19,9; 4,6; 6,7; 7,8; 10,10; 9,12; 18,15 

18 

Hydrabad 

10,2; 9,4; 18,9; 10,6; 9,8; 18,9; 4,10; 6,11; 7,12 

4 

Lucknow 

2,1; 9.4; 8.9; 10,6; 12,9; 2.10; 9,12; 8,13 

9 

Bangalore 

4.1; 6,2; 7.3; 8.6; 9.8; 10,10; 10,11; 10,12; 19,13 

6 

Jamshed. 

8.2; 9,4; 10,6; 10,7; 10,8; 19,9; 8.10; 9,12 

1 

Delhi 

10,1; 12,4; 2,6; 9,8; 8,9; 10,10; 12,13 

12 

Bombay 

10.1; 10,2; 10,4; 19,9; 9,6; 10,7; 19,9; 9,10; 



10,11 ; 19,13 
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5*2*3 Da'ta about, t-ho demand requirements 

Table 5*3 Demand at various nodes 


Demand node 

Demand, time (day) 

19 

Bhopal 

10.?; 1?.6; 10.9; ?.10; 4.11; 8.15; 7.14 

15 

Ahmadabad 

2.6; ?.7; 10.9; ?.ll ; 8.15; 7.1? 

7 

V i jaywada 

20.6; 10.8; ?.9; 1?.12; 1?.15; 10.14 

8 

Madras 

1?.7; 10.8; 20.9; 10.11; 10.15; 1?.14 

14 

Manmad 

?.4; 1?.8; 2,9; ?.10; 2,15; ?.14 

16 

Nagpur 

8.?; 7.6; 8,9; 7.10; ?.15; 10,16 

? 

Calcutta 

10.4; ?.6; ?.9; 10.12; ?.15; 1?,16 

20 

Indore 

?,4; 10.7; 1?,9; 10,10; 10.15; ?.14; 4,16 

5 

Kanpur 

10,4; ?,?; 4,6; 10,9; 1?.10; 10,15; ?.l? 


The empty wagons made available in the entire network in 
the second time period is varied for both problems. Solution is 
obtained for four cases of each problem. 

1. 17? 2. 200 J. 22? 4. 2?0 

5* 3 OUTPUT 

The output solution is volumnious. We get following 
different outputs. 

1. Distribution of number of trains delayed vs. amount of delay 
Percentage frequency of delays vs. amount of delays. 
Percentage of delay vs. amount of delays 


Average delay. 
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2. Routes followed by the trains and schedule of the trains 
separately for both goods movement and empty wagon movement. 
It IS expressed as nodes on the path and the starting time 
from each node on the path which is stated taking intc 
consideration idle time if any at a node. 

5. Load on each arc, each day. 

The maximum load on the arc. and the time at which it occurs 

4. The cost of scheduling for both goods trains and empty trams 
in each time period. And total cost of the schedule 
generated . 

5. 4 FINAL RESULTS 


Table 5.4 
Problem 1 


Number of 

Total Cost 

Average 

Empty Wagons 

of Schedul i ng 

Delays 

17? 

5?.886 

?.l? 

200 

??.615 

?.2? 

22? 

5?, 096 

4.90 

2?0 

J5.957 

4.?5 




4 ? 


Table 5* 8 
Probl em2 


Number of 

Total Cost 

Average 

Empty Uagons 

of Schedu 1 i ng 

De lays 

175 

38.786 

3.49 

200 

38.416 

5.37 

225 

38.020 

3.30 

250 

36.706 

5.01 


This results are plotted in the figure 5.2 and figure 5.3 









Totftl 
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S* 5 DISCUSSION OF THE RESULTS OBTAINED 

Fig. . 2. gives a graph of total cost of transportation 
vs. number of empty wagons available. One can very easily notice 
that the total cost of the transportation decreases with the 
increase in the number of empty wagons available. It is simple to 
interpret it. When the number of empty wagons available is 
inadequate, then these need fetched from far off places before 
these can be used for goods transportation, which increases the 
total cost of the transportation. 

Figure 5.3 a graph of Delays in the movement of the 
wagons. In this graph also it can observed that there is decrease 
in delays with increase in the availability of the wagons. Very 
similar explanation seems to be valid for this behaviour also, 
i.e. more the number of empty wagons available, lesser will be the 
need to bring them from over a long distance. However a tradeoff 
IS possible for the decision maker between higher cost of more 
trains and the lesser delays and transportation costs. 

If we compare the results obtained for two problems then 
we find that the problem 2 in which we incur more total cost of 
transportation has correspondingly lesser amount of delays taking 
place and vice versa. However in problem 1 we find that as the 
number of empty wagons available is decreased' from 200 to 175 
there is a slight decrease in the delays (even though the cost of 
transportation has increased). This slight discrepancy probably 



due to the fact that our 


"single pass" 


n nature. 


heur I s 1 1 c is 



CHAPTER 6 
CONCLUSION 
AND 

DIRECTIONS FOR FUTURE WORK 


6. 1 CONCLUSIONS 

The computational experience with the algorithm that has 
been developed seems to be reasonably well. The moderately sized 
problem, has been solved on HP 9000, in approximately 30 secs of 
the CPU time. 

We feel that the graphs that are obtained, can be very 
useful for Railways for taking decisions regarding the investments 
to be made in the critical assets like wagons. In reality usually 
there is some costs associated with the delays e.g. Overtime paid 
to the employees, more fuel and other operating costs for engines. 
If the cost of delays can be quantified, then it should be 
possible to find the sum of cost of transportation and cost of 
delays. It will help Railways for getting a idea of approximate 
investments desired so as to have the minimum net cost of trans- 


portat ion. 
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6*2 DIRECTIONS FOR FUTURE WORK 

After working on the problem of Railway Wagon Scheduling 
we have realized various dimensions of this area. There is 
scope for continuing work on this problem area. 

It will be useful to work in this area to develop a 
methodology for dealing with multicommodity problem. It may be 
enough to deal with few goods. This is because 80% of the total 
freight movennent is due to few types of goods like steel. food 
grain, fertilizers, cement. Attempt could be made by modifying 
methodology so as to incorporate multicommodity min-cost flow and 
multidimensional assignment problem. 

Given the difficult nature of the problem it is felt 
that heuristics need to be developed. Also heuristic must be of 
the "multipass" type. So that each time we are able to improve on 
the solution of the previous iterations. 

The scope of the problem can be widened so as to a 
decision support system for other important investment areas of 
Railways e.g. Rails of the Railway lines. Certain rail routes are 
heavily loaded, where as others are scarcely utilized. The study 
of this aspect will be useful in taking decisions regarding 
increasing the capacity of the routes. 

It will be desirable to take a real life data. Also, an 
attempt needs to be made to solve at least one or two problems 
optimally. It will give us the yardstick to compare the 


perform- 
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ance of the heuristic approaches developed for such problem. 
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APPENDIX A 


Here we give a formulation for pick and delivery problem 
with time windows from Solomon C?1D. Then this model is used to 
develop a taxonomy of routing and scheduling problems with time 
windows . 


Let there be n customers indexed by i . Associate to the 
pickup location of customer i a node i and to his delivery 
location a node n+i. Also associate to the depot node 0 and 2n + 
1. This creates a clear distinction between customers. their 
associated locations and the nodes of the network. Note, however, 
that different nodes may correspond to the same physical location. 
Therefore. N = C1.2....n, n+1 . n+2 . . . . 2n+l > is the node set for our 
network, and P = P^U P . P^= Cl,2,..,n>, P = Cn+1 , n+2,..,2n} is 

the set of nodes other than the depot nodes. We will index both N 
and P by u and w. Customer i demands that d^ units be shipped from 
node i to node n+i . Next, let Ca^.b^] denote the pickup tie 
window for customer i and let Ca ^ denote his delivery time 

window. Let V = £1.2 IVI) be the set of vehicles to be routed 

and scheduled; index this set by v. Let D be the capacity of each 
vehicle. Let also Ca .b ] denote the vehicles’ dparture time 

O O 

window from the depot and their time window for 

arrival back at the depot. 
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For each distinct u,w in N, let t and c represent the 

uv uv 

travel time and the travel cost form u to w, respectively. Retain 
only those arcs <u,w) for which (1) a + s + t ^ b , where s 

U U UV V u 

IS the service time at customer u, and either (2> d + d < D. u.w 

U V 

e P , u w. or (2') d + d ^ D, u.w e P , u w. Finally. 

u — r» V- n 

let K be the fixed cost associated with each vehicle, incurred if 
this vehicle is utilized. 

Three types of variables are used in the mathematical 


formulation : binary flow variables , u. w e N. v e V. u ^ w, 
time variables T . . T . u e P. v e V. and load variables Y . 

u o 2n4.i 

u e P- Let the binary decision variable equal to one if 

vehicle v travels from node u to node w. and equal to zero 
otherwise, v e V. and u. w e N, u w. Let Tu be the time at 
which service at customer u begins, u € P. Let also be the 

time at which vehicle v leaves the depot and be the time at 

which it returns to the depot, v € V. Next, let Y^ be the total 
load on the vehicle just after it leaves node u. u e P. It is 
assumed that the vehicles depart empty from the depot. i.e.. 


0 


constraints . 


These are decision variables with nonnegativity 


Finally, let f <T ) denote the penalty associated with 

beginning service at node u at time T^ , u € P, and g < 2 n+i o 

* ^ o T>rkii'i“p of the indicated duration 

be the penalty associated with a rou 

. ^ ^ V These penalties must be 

when traversed by vehicle v ^ 

r-riQtA c . We are now in a 

expressed in the same units as t 
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position to present the pickup and delivery problem with time 
windows formally. The mathematical formulation is : 


Min v^ieN 5>eN + BiSP fu(Tu) + X>ev ) <1> 




^8 


X'' =1=>Y - d =Y. 

uv u w V 


> 

X 

= 1 ^ Y + 

d 

Ov 

o 

V 

Y 

o 

= 0, 0 < Y 

< D 

x'' 

uw 

b 1 nary , 



u e P , w e P , 

V € V (1^> 

weP*, veV (16) 

u € P^ (17) 

u.weN, veV (18) 


We seek to minimize the sum of the total cost, the total 
penalty associated with servicing customers too early or too late, 
and the total penalty associated with routes exceeding a given 
duration. Constraints (2)-{5> and <18> forma multicommodity 
min-cost flow problem. Constraints <6) ensure that the same 
vehicle, v, visits both u and n+u. Constraints (7) are precedence 
constraints which force node u to be visited befor node n+u. 
Next, constraints {8>-(10) describe the compatibility requirements 
between routes and schedules, while constraints are the 

time window constraints. Finally, constraints (14>-(16) express 
the compatibility requirements between routes and vehicle loads, 
while constraints (17) are the capacity constraints. 


Note that the formulation includes a route duration 
restriction (b - a >. Note also that the formulation does not 

2 n+ 1 O 

include subtour elimination. 

L.et us now examin how models for different problem 
variants can be obtained from this formulation. The above 
forulation is the well known mu Itivehicle, man y to many 
dial-a-ride problem with time windows. If we let IVI = 1 and 
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elirrtin8.te constraint (6) we obtain the single vehicle variant. 

Redefine now N = P U CO.n+U. P = p'*' = p"= C1.2....n>. 

To obtain the VRPTW we delete constraints {6>. <7) and (15). The 
constraints (8), as well as (14) ensure subtour elimination. If 
the VRPTW involves only pickup, define d > 0. For delivery 

V 

VRPTW, d ^ 0, and in constraint (17), set Yo = D and eliminate 

w 

Y < D. The multiperiod VRP can be viewed as VRPTW where the 

U 

customers are replaced with request for service and for each 
request there is a time window. 

Furthermore, different schortest path problems with time 
window can be obtained from the above formulation. For example, 
constraints (1), (5)-(5), (8)-(15), (18) and IVI = 1 constitue one 

shortest path problems with capacity, precedence and time window, 
can be obtained similarly. 
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PROGRAM FOR RAILWAY WAGON SCHEOULING 


PROGRAM heurilO (INPUT. OUTPUT); 

C Trying to determine costs. Used for solving problems 
with restricted capacity of arcs, so as to make mincost infiJ 
LABEL 25 . 50 , 7? , 100 , 200 , 2^0, 500 . 550.575,580.550.575.580,600; 

CONST mmax=100; nmax=200; nodmax=75: arcmax=150; 

T_hormax = 120; C Max Time horizon (unit time) 1 

TYPE 


C mainly used for the proce simplex ) 
matrix = ARRAY Cl..mmax] OF REAL; 
matrint= ARRAY Cl..mmax] OF INTEGER; 
row_matrix = ARRAY Cl..nmax] OF REAL; 
opt imi sat 1 ontype = ( maximisation, minimisation); 
bigmatrix = ARRAY Cl . .mmax , 1 . . nmax] OF REAL; 
inequal array = ARRAY [1.. mmax 3 OF CHAR; 

C mainly used for defining node’s data > 
nodetype = 0..nodmax; 
nodeset = SET OF nodetype; 
nodrec = RECORD 

dema_supi : ARRAY [1 . . T_hormax] OF REAL; 
cum_dem_sup_g ; REAL; 

requ_avai : ARRAY CO . . T_hormax3 OF REAL; 
cum_dem_sup_w : REAL; 

END; 

nodaray = ARRAY Cl..nodmax] OF nodre^c; 

C mainly used for the defining arc's data > 
arc_rec = RECORD 

nodel , node2 : nodetype; 
dist,cost_g. cost_w. 
capa, tmp_capa : REAL; 
deleted. enter : BOOLEAN; 

res_cap,load : ARRAY Cl - .T_hormax3 OF REAL; 
t im_s t , t 1 m_dt : REAL; 

END; 

arcaray = ARRAY Cl..arcmax] OF arc_rec; 



£ used for procedure f i nd_ass i gnment } 
dem_nod = RECORD 

s rtodelypej 
dem : REAL; 
t im ; 1 . . T_hormax ; 

END; 

ass_dem_nd = ARRAY [l..nodmax] OF dem_nod; 

8up_nod = f^LCOf^U 

nd : nodetype; 
sup : REAL; 
t im ; 1 . . T_hormax ; 

END; 

a88_sup_nd = ARRAY Cl,.nodmax3 OF sup_nod; 
sl._ass_typ = ARRAY C 1 . . nodmax , 1 . . nodmax3 OF REAL; 
arcset = SET OF INTEGER; 
arcarr = ARRAY Cl..arcmax3 OF INTEGER; 
del_type = ARRAY [0 . . T_hormax3 OF INTEGER; 
wtd^de lay_typ = ARRAY [ 0 . . T_hormax3 OF REAL; 

VAR C mainly used for simplex procedure > 
mm, mml , nnl , n 1 , 1 i , j j , kk : INTEGER; 

Cm ; r ow^^ma t r i x ; 

Am : b 1 gma t r i x ; 

inequal itym : i nequalarray ; 

bmi , bmo : matrix; 

ans : CHAR; 

pariod_cost : real; 

inp,OUt,OUt^ : TEXT; 

C u^ed in mincost network > 

tt,T_, per iod_lenth,T_hor : O..TJiormax; 
per iod_no : 1 . , 10 ; 

no^nod , no^arc, nod , nodi , nod2 : nodetype; 
arc 2 arcaray; 
node : nodaray; 

degree: ARRAY Cl-.nodmax3 OF INTEGER; 
ds : REAL; 

mi ncos t_f eas i bl e : BOOLEAN; 

C Note that since all the arcs undirected i.e. two way 
indegree = outdegree = degree /2 > 
cumsum : INTEGER; 

point, fmovptr. rmovptr : ARRAY Cl..nodmax3 OF INTEGER; 



f trace, rtrace : ARRAY Cl..arcmax3 OF INTEGER; 

C used for procedure f i nd_as 8 i gnment } 

dnn.snn. ino.jno, s_nod . cur_nod . d_nod : nodetype; 

dmn : ass_dem_nd: spn : ass_sup nd; 

slak. assign i sl_ass_typ; 

dem_nod_set : nodeset; 

arc_flo_set : arcset; 

cn , f t , arc_s tor : INTEGER; 

assl ,ass 2 .ass , mi n_arc_capa , prev_mi n_arc_cap . flo_sent. 
aval l_t im. mot 1 on_t im, travl_t im. s lack_t im. de lay . rk REAL; 
s_tim.d_tim. time O..T_hormax; 

ddeg.sdeg. spt.dpt ; ARRAY Cl..nodmax3 OF INTEGER; 
arcs_aray : arcarr; 

C Empty wagon scheduling } 

no. i store . j store . kstore . ptr : INTEGER; 

excl_arc : arcset; 

max load : REAL; passed, tried : BOOLEAN; 

C output information > 
delay_freq : del_type; 
wtd_delay : wtd_de lay_typ ; 
tot_de lay . total_cost : REAL; 
out! , out2 , out5 , out4 : TEXT; 

FUNCTION smaller <a,b ; REAL) : REAL; 

BEGIN 

IF a < b 

THEN smaller ;= a 
ELSE smaller := b; 

END: 

PROCEDURE simplex 



(m.n: INTEGER: C: row_matrix: A.-bigmatrix; inequal ity: inequal array : 
VAR b_oimatrix; VAR dec i_vars : arcset ; VAR feasibility : BOOLEAN) : 

LABEL ^^,10,20; 

TYPE set_^of^var - SET OF l..nmax; 

VAR optimisation : opt imi sat lontype ; 

set_of_slack_var. set_of_art i f i c ial_var . set_of _al trn_optm_ 

set^of var 5 


b: matrix 


soln t 





reducedcost : row_matr i x ; 

b_bar, A_onter i ng . CB . PI . rat i o i matrix; 

basic_var : array [l..mmax] of INTEGER; 

B_inv : ARRAY C 1 . .mmax . 1 . .mmax] OF REAL; 

non_ba8_var; ARRAY tl..nmax] OF INTEGER; C only ml..nt used 3 

1 .ml ,nl ,nt, iterati on no , new_n , o 1 d n , deleted, print no.j, count , 

pivot_col . pc , pivot_row, basic_col. 

no_of_hyphonsl , no_of_hyphons2 : INTEGER; 

min_ratio, mi n_reducedcost ,max_reducedcost . Z , product : REAL; 

PROCEDURE Find_max_reducedcost(VAR max : REAL; VAR indexcolumn : INTEGER) 
VAR il : INTEGER; 

BEGIN 

max := reducedcostCm+l 3 ; indexcolumn := m+1 ; 

FOR 1 1 : = m+2 TO new_n DO 
IF <max <= reducedcost C 1 1 3 > 

AND not (non_bas_varC i 1 3 in set_of_altrn_optm_soln) 

THEN BEGIN 

max := reducedcost C 1 1 3 ; 
indexcolumn := il 

END; 

END; 

PROCEDURE Find_min_reducedcost<VAR min : REAL; VAR indexcolumn : INTEGER) 
VAR il : INTEGER; 

BEGIN 

min r educedcoB t Cm+1 3 ; indexcolumn := m+1; 

•FOR il .•* m+2 TO new_n DO 
IF (min >= reducedcostC 1 1 3 ) AND 
not<non_bas_varC i 1 3 IN set_of_al trn_optm_soln) 

THEN BEGIN 

min := reducedcostC i 1 3 ; 
indexcolumn := ili 

END 

END; 

PROCEDURE Find_min_ratio(VAR min : REAL; VAR indexrow ; INTEGER ); 

VAR indexro : ARRAY Cl..mmax3 OF INTEGER; 

C Stores all the ratios which are minimum & also equal 3 

stackptr*il : INTEGER; 



clegeneracy_f ound : BOOLEAN; 

PROCEDURE R®so 1 vB_degeneracy ; 

C j 

VAR column. i2 : INTEGER; 

tempratio, minratio : REAL; 
dogeneracy_resolved : BOOLEAN; 

BEGIN C Resolve degeneracy } 

degeneracy_resolved := FALSE; column := 0; 
WHILE ( NOT degeneracy_resol ved ) AND ( column <= m > DO 
BEGIN 

column := column + 1; 

minratio := 10000; C initialise } 

FOR i2 := 1 TO stackptr DO 
BEGIN 

IF A_enter 1 ngC 1 ndexroC 1 2 ] ] > 0 

THEN tempratio := B_i nvC i ndexrot i 2] , column] / 

A_enteringC indexroC i23] 

ELSE IF A_enteringC indexroC i2]3 = 0 
THEN tempratio ;= 10000 
ELSE tempratio := -10000; 

IF minratio > tempratio 
THEN BEGIN 

minratio := tempratio; 
indexrow := indexroCi23; 

END; 

END; 

degeneracy resolved := TRUE; 

END C while 1 

END; C Resolve degeneracy 1 

BEGIN C Find min_ratio > 

dageneracy_f ound :* FALSE; 
mi n ! = -1 ; 

FOR i 1 : = 1 TO m DO 
IF ratioCil] >= 0 

THEN IF (min < 0) OR (mm > ratioCil]> 

THEN BEGIN 

rnin := ratioCills 
i ndexrow : = i 1 ; 
degeneracy_f ound := FALSE; 



stackptr : = i . 



indexroCstackptr] •= ii- 

£NU 

ELSE IF min = ratioCil] 

THEN BEGIN 

degeneracy_found := TRUEj 
stackptr := stackptr + 1; 
indexroCstackptr] : = ilj 

END; 

IF degeneracy_found = TRUE 
THEN Reso 1 ve_degeneracy ; 

END: C Find min ratio } 

PROCEDURE Interchange_suff ixes_of_variables<suff 1 x 1 .suff 1 x 2 : INTEGER); 
VAR temp : INTEGER; 

BEGIN 

temp : = non_ba3_varCsuf f 1 xl ] ; 

IF suffix2 <= m 
THEN BEGIN 

non_bas_varCsuff 1x1 ] := bas i c_var[ suf f i x2] ; 
bas 1 c_var[ suf f 1 x2] := temp 

END 

ELSE BEGIN 

non_bas_varCsuf f i xl ] ;= non_bas_varCsuf f ix2] ; 
non_baa_var[suff 1 x 2 ] := temp; 

END; 

END: 

PROCEDURE Interchange(VAR x.y ; REAL); 

VAR temp ; REAL; 

BEGIN temp := x; 

X : = y ; 

y : = temp 

END; 

PROCEDURE Interchange_columns_of_var lables 
(enter ing_bas i 8 , 1 eav ing_bas 1 s ; INTEGER); 

VAR il . INTEGER; 



VAR temp : REAL: 
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BEGIN 

FOR 1 1 := 1 TO m DO 
BEGIN 

temp := AC i 1 , entering_bas is] ; 

ACil .entering_basis] := A[ i 1 . 1 eav i ng_bas i s ] ; 
AC i 1 , ieaving_basis] .-= temp 

END 

END; 


BEGIN £ Procedure simplex } 

REWRITE £ out , ' outpro ) e . dat * > ; 
pr 1 nt_no : = 1 ; 
ml : = m + 1 : 

C To find out total no variables (nt) 1 
nt := m+n; 

FOR 1 : = 1 TO m DO 
IF 1 nequal i ty C 1 ] = ’>' 

THEN nt := nt + 1; 
new_n : = nt ; 

optimisation : = minimisation; 

C Columns indexed m+1 to new_n are for non_bas_var 

and columns indexed 1 to m are basic_var 1 
C Initialise non_bas_var and ba8ic_var 1 

• FOR i;= 1 TO nt-m DO non_bas_varCm+ 1 ] := i; 

FOR 11 = 1 TO m DO basic_varCi3 ;= nt-m+i : 

CASE optimisation OF 

maximisation i min_reducedcost :=-l ; 

minimisation ; max_reducedcost i= 1; 

END; 

n 1 i=m+n: 

set_of_siack_var := C3; set_of_artif icial_var i= C]; 
set__of__al trn_optm_soln ;= C]; 

C Read A & b > 

count := 0; 

FOR 11=1 TO m DO 
BEGIN 

IF (nl+1 <= nt) 
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THEN FOR nlFl TO nt DO 

AC 1 , J 3 0; 

CASE inequal i tyC I ] OF 
*<’ : BEGIN 


+ [bas i c_varC i 3 3 


+ Cnon_bas_varCnl+count3 3 ; 


•f Cbas 1 c varC i 3 3 


CCi3 :=0; 

set^of^s lack^^var := set of slack var 


END; 

: BEGIN 

count : = count -f 1 ; 

AC 1 ,nl+count3 := -1; 

CCnl+count3 0; 

set_of_s lack_var r = set_of_s iack_vai 

CASE optimisation OF 

maximisation : CCi3 := > -1000; 

minimisation : CCi3 := £ M 3 1000 
END; 

set of artificial var := set of artj 


END.; 

*=* : BEGIN 

CASE optimisation OF 

maximisation : CC i 3 := £~M >“100( 
minimisation : CC 1 3 := £ M } 1001 


END; 

set__of_ar 1 1 f 1 c i al_var : = 

80t_of_^art I f i c lal^var 4- Cbas i c^varC i 3 3 

END; 

END; £ case inequal i tyC i 3 > 

CBCi3 := CCi3; 

END; 


C Initialise coefficients of basic variables & B_inverse as unit matrix 

FOR 1 : = 1 TO m DO 
FOR j * = 1 TO m DO 
IF i-i 
THEN BEGIN 

AC i . j 3 != 1 ; 

B invC i , j ] : = 1 
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END 

ELSE Or.lilN 

At 1 , j 3 := Oj 
0_invC 1 . j 3 := 0 

END; 

C Initially find the SIMPLEX Multipliers } 

FOR j := 1 TO m DO 

BEGIN 

PIC j3 := 0; 

FOR 1 :=1 TO m DO 

PICj3 := PIC)3 + CBCi3 * B_invCi.j3 

END: 

C Initially Find Reducedcosts for all variables 3 
FOR j : = 1 TO new_n DO 
BEGIN 

product := 0; 

FOR i : = 1 TO m DO 

product := product + PICi3 * ACi,j3; 
reducedcost C J 3 := C[j3 “ product 


END; 

£ Initially take b_bar = b } 

FOR 1 := 1 TO m DO 

b_barCi3 := bCi3; 
iteration_no r= 0; deleted 0: 

CASE optimisation OF 

maximisation : F i nd_max_reducedcost ( max_reducedcost 


pi vot_col > J 


minimisat ion 


Find min_reducedcost( mi n_reducedco8t 


pivot_col ) ; 

END; 


5 1 REPEAT 

iteration no ;= iteration_no + 1 ; 

C Compute the updated coeffit.ents of entering column ) 

FOR is® 1 TO m DO 

BEGIN 

A_enter ingC i 3 := Oj 

FOR j := 1 TO m DO ' 

A_enteringCi3 := A_enter ingC i 3 + B_.nvCi.j3 

* AC j ,pivot__col 3 


END; 



<= 0.001 


C Apply minimum ratio rule for deteriming leaving variable 
F ind rat los } 

FOR 1 :=1 TO m DO 
BEGIN 

IF abs < A_enter 1 ngC 1 3 - ro’und(A_enter ingC i 3 > > 


THEN A_enter 1 ngC i 3 := round( A_ent er i ngC i 3 > ; 
IF A_enter ingC i 3 > 0 

THEN ratioCj] := b_barCi3 / A_enter i ngC i 3 
ELSE IF A_enter ingC 1 3 = 0 

THEN ratioCi] := 10000 
ELSE ratioCil := -10000; 

END: 

F i ncl_fr)i n_rat 1 o{ min_ratio. pivot_row >; 

IF min_ratio < 0 
THEN BEGIN 

WRITELN (out. 'This problem has 


unbounded solution ’>; 


WRITE ('Unbounded solution !!. 


Program aborted . . ’ ) ; 

GOTO 20 


C Update B 
FOR 


END: 

inverse matrix 3 

jjal TO m DO (Update pivot row firstl 
IF A enteringCpivot_row3 > 0 

THEn"b invCpivot_row.)3 := B_i nvCp i vot_row . j 3 / 


A_enteringCpivot_row] ; 

FOR i s = 1 TO m DO 


{ Update all other rows 3 


IF 1 <> pivot_row 
THEN FOR J TO m DO 

B invC i . j 3 s = B_invC i . ) 3 


- 6 invCpivot_row, j 3 


* A_enter ingC i 3 ; 

{ Find b_bar 3 

FOR i s =1 TO m DO 
BEGIN 

b_bar C i 3 s ~ 0 • 

FOR }!~ 3 TO m DO 
b barCi3 b_barCi3 


+ B invCi,j3 * bCi3 


END; 

basic col s- pivot_rows 



£ Find SIMPLEX multipliers 'Pi's } 

COCbas 1 c_col ] := CCpi vot_col ] ; 

r OR J : = 1 ro m DO 

BEGIN 

PICj] := 0; 
for 1 :=1 TO m DO 

PICj] := PICj] + CBCi] ♦ B_invCi.j3 

END; 

C Find Reducedcosts for all variables. 

note that reducedcost will automatically be = 0 

for basic_var } 

FOR j := 1 TO new_n DO 
BEGIN 

product := 0; 

FOR 1 : = 1 TO m DO 

product := product + PI[i] * ACi,)]; 
reducedcostC j 3 := C[ j 3 - product 

END; 

Int erchange_suf f i xes_of_var i abl es (p i vot_col , bas i c_col ) ; 
Interchange<CCpi vot_col ] .CCbas i c_col ] ) ; 

I nterchange_co lumns_of_var lables ( pi vot_col , bas i c_col ) ; 
Interchange < reducedcost Cpi vot_co 1 3 , reducedcost C bas i c_col ] > 
pc := non_bas_varCpi vot_col 3 ; 

IF non_bas_varCpi vot_col 3 IN set_of_art i f i c i al_var 
THEN BEGIN 

£ Delete the artificial variable from subsequent considerations 
once it has left the basis > 

old_n :* nt - deleted; 

I n t e rchange_su f f i xes_of_var i ab 1 es ( p i vo t_co 1 , o 1 d_n ) ; 

I nter change <C[pivot_col 3 .CCold_n3 > ; 

Interchange_columns_of_var iables(pi vot_col , old_n) ; 

Int e r change (reducedcostCpi vo t co 13. reducedcost Col d_n 3 ) ; 

deleted deleted + 1: 

END; 

feasibility := TRUE; 

FOR i:= 1 TO m DO IF b_barCi3 < 0 
THEN BEGIN 

WRITELN (out, 'Now problem is primal infeasible, 

because one rhs constant b_bar is -ve'); 

CASE optimisation OF 



/ 1 

maximisation : FOR j:= ml TO new_n DO 

IF not()'IN set_of_s 1 ack_var ) 

and not<j IN aQL_of_.Tirtif icial_var> and 

(reducedcostC j ] > 0> 

THEN BEGIN 

WRITELN 

(out, 'It IS also dual infeasible. Hence can’'t be solved further'); 

WRITE 

('Infeasible solution !!. program aborted 

feasibility := FALSE; 
GOTO £0; 

END; 

minimisation : FOR ):= ml TO new_n DO 

IF not(j IN set_of_slack_var > 

AND not(j IN se t_of_art i f i c i a l_var ) and 

(reducedcostC j 3 < 0) 

THEN BEGIN 

WRITELN 

(out. 'It IS also dual infeasible. Hence can''t be solved further'); 

WRITE 

('Infeasibible solution 11, program aborted . . ' ) ; 

feasibility ;= FALSE; 
GOTO ZQ 
END; 


END; 

WRITELN (out , 'However you may solve it further, 
by dual simlex method'); 

WRITE ('Primal infeasible !!, program aborted ..') 
GOTO 20; 


END: 

new_n ;= nt - deleted; 
CASE optimisation OF 

maximisation : BEGIN 


(max__reducedcost , pi vot_col ) ; 
(max_reducedcost) ) <= 0.001 
round(max reducedcost ) ; 


Find_max_reducedcost 
IF abs( max_reducedcost 
THEN max_reducedcost := 
IF max reducedcost = 0 


round 



then set_of_altrn optnri_soln : = 

set_of_altrn_optm_3oin + Cpc] CLnon_bas_var[pivot_col]]}; 

END; 

minimisation : BEGIN 


<mi n_reducedcost ,pivot_col ) ; 
roundCmi n_reducedcost ) > <= 0.001 

round(min reducedcost > ; 


F 1 rid_mi n_reducedcost 
IF abs( mi n_reducedcost - 
THEN min reducedcost := 


IF mi n_reducedcos t = 0 
THEN set_of_aItrn_optm_soln : 

= set_of_al trn_optm_soln + [pc] C Cnon_bas_varCp i vot_co 1 ] ] } ; 

END 

END; 

UNTIL ( max_reducedcost <= 0 > OR ( min_reducedcost >= 0 ) ; 

C Furtfier minimisation of objective function is not possible 
i.e. when optimum solution is obtained }; 

IF set_of_aI trn_optm_soln <> C] 

THEN BEGIN 

WRITE <out,'Set of alternate optimum soln contains ; ’); 

FOR i := 1 TO nt DO 

IF 1 IN set_of_altrn_optm_soln 

THEN WRITE (out . ’X’ , i :2) ; WRITELN (out); 

END; 

10: C Compute the value of minimum value of objective function! 

Z := 0; 

, FOR i ! = 1 TO m DO 

Z := Z + CBCi] * b_barCi]: 

C Print output } 

CASE optimisation OF 

maximisation : WRITELN(out . 'Maximum objective function value 
is Z = ',Z:6:2>; 

minimisation ; WRITELNCout . 'Minimum objective function value 
is Z = ' ,Z;6;2> 

END; 

FOR i := 1 TO m DO 

b_oCi3 := 0; C initialise before returning output J 
WRITELN<out , 'Opt imum solution is '>: 

FOR i ;= 1 TO m DO 



BEGIN 


IF NOT ( <bas i c_varC 1 ] IN set_of_s lack_var ) OR 
(bas 1 c_varC i 3 IN se t_of_ar t i f i c i al_var ) ) AND <b_bar[i] > 0.001) 

THEN BEGIN 

WRITE (out.* X’ .basic_varCi]:2. ' = * . b_barC i 3 : 4 : 2 > ; 
deci_vars := deci_vars + [bas i c_varC i ] 3 : 
b_oC bas i c_var[ 1 3 3 := b_barCi3: 

END; 

END; WRITELN<out) ; 

FOR I := 1 TO m DO 

IF ( bas 1 c_var C 1 3 IN se t_of_art i f i c i al_var ) AND (b_barCi3 <> 0 ) 
THEN BEGIN 

WRITELN (out. 'Since artificial variable is present 

in the optimum solution *); 

WRITELN (out, 'PROBLEM IS INFISIBLE'); 

WRnELN(out.* X* .basic_varCi 3 :2 , ' = ’ . b_barC i 3 : 4 : 2 > ; 
WRITELN ( ' Inf 19 ibi 1 1 ty !!. program aborted 
feasibility := FALSE; 

GOTO 20; 

END 

ELSE feasibility := TRUE; 

CASE optimisation OF 

maximisation .• IF (max_reducedcost = 0) 

AND not(non_bas_varCpivot_col 3 IN set_of_altrn_optm_soln) 

~ THEN BEGIN WRITELN (out. 

'Finding out another alternate optimum soln ');GOTO ^ END; 

minimisation : IF (min_reducedcost = 0) 

AND not ( non_bas___varCpi vot_col 3 IN sat_of_altrn_optm_soln) 

THEN BEGIN WRITELN (out . 'Finding out 

another alternate optimum soln ');GOTO 5 END, 

END; 

WRITE ('Complete solution has been written in the 

file outproje.dat’); 

20: READLNs {CLOSE ( out ); > 

END; C Procedure Simplex 3 

FUNCTION shortest_dist ( s. t : nodetype): REAL: ^ 



VAR maxdist, temp_dist : REAL; 

Naet, temp, perm : nodeset; 
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shordist : ARRAY [l..nodmax] OF REAL; 
j.j, arc_no : INTEGER; 
stoi'e, stoi'l : nodetype; 

function min_temp_dist_nod : nodetype ; 

VAR minimum ; REAL; k : nodetype; 

BEGIN 

minimum := maxdist; 

FOR k := 1 TO no_nod DO 
IF NOT <k IN perm) 

THEN IF shordistCk] < minimum 
THEN BEGIN 

minimum := shordistCk]; 
min_temp_di st_nod := k; 

END; 

END; C mi n_temp_d i s t_nod > 

BEGIN C function shortest dist } 

Nset := Cl..no_nod]; 

perm ;= Cs3; temp := Nset - Cs]; 

IF NOT (<B IN Nset) AND (t IN Nset)) 

THEN BEGIN 

WRITELN ('Node out of network data range !!'); 

GOTO 500; 

END; 

£ Sat shortest dist for other nodes = sum of all dist.s in network > 
maxdist := 0; 

FOR ii= 1 TO 2*no_arc DO 

maxdist := maxdist + arc[i].dist: 

FOR i != 1 TO no_nod DO 

shordistCi] := maxdist; 
shordistCs] .*= 0; 
store := s; 

WHILE NOT <t IN perm) 00 
BEGIN 

FOR pointCstore] TO pointCstore+1 3 -1 DO 

BEGIN C Get the node in the adjuscency list of node 'store' 
called 'storl', update its shordist if required > 
arc_no := ftraceCj]; 

temp_dis.t := shordi stCstore] + arcCarc_no3 . di st ; 



storl := arc[arc_no] . nodeZ ; 

IF storl IN temp 

1HFN BhordistCstorl 1 i« omallor ( shord i s t [storl 3 » 

temped ist) ; 

END; C for adjuscent nodes update shordist } 
store := mi n_temp_di st_nod ; 
temp := temp - [store]; 
perm := perm + [store]; 

END; C while not (t in perm) DO } 
shortest_di st := shordistCt]; 

END; C function Shortest dist > 

procedure Find_Assignment < nodn.nosn : nodetype; 

} dn ! ass_dem_nd; sn : ass_sup_nd; 

slack : sl_ass_typ; 

VAR assignment : sl_ass_typ>; 

VAR ip. iP. ro.co, count ; INTEGER; 

assigned : ARRAY [ 1 . . nodmax , 1 . . no'dmaxD OF BOOLEAN; 

totslack. small : REAL; 


procedure Max__slack <VAR row, col: INTEGER) : 

VAR i.j ; INTEGER; max : REAL; 

PROCEDURE braketie; 

VAR penaltyl, penaltyZ : REAL; 

FUNCTION next_max_in_row (iro : INTEGER/ : REAL, 

VAR k : INTEGER; nextmax : REAL; 

BEGIN 

nextmax := 0; 

FOR k 1 = 1 TO nodn DO 

IF <na><tmax < alackC. ro . K J ) AND (alackC. ro . k] < aax> 

AND NOT ass i gnedC 1 ro , k ] 

THEN nextmax ;= s lackC i ro , k] ; 
next max_in_row := nextmax; 


END; 

FUNCTION next_max_in_col <ico 
VAR k ; INTEGER; nextmax 

BEGIN 


INTEGER) : REAL; 
REAL; 


nextmax := 0; 


FOR k := 1 TO nosn DO 
IF (nextmax < s lackCk . i co3 ) AND 


(si ack[k 


co] < max) 


/ 
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and not ass I gnedCk , 1 CO] 

THEN nextmax := slackCk , i co] ; 
naxt_max_tn_col := nextmax i 

END; 

begin C brakatie } 

IF j = col 
THEN BEGIN 

penal tyl .- = max - next_max_i n_row( i ) ; 
penal ty2 := max - next_max_i n_row( row) ; 

IF penalty! > penaltyZ 
THEN row : = i ; 

END; 

IF i = row 
THEN BEGIN 

penalty! := max - next_max_i n_co 1 ( ) ) ; 
penaltyZ := max - next_max_i n_co 1 < co! > ; 

IF penalty! > penaltyZ 
THEN col := j ; 

END; 

END; C braketie } 

BEGIN C procedure max_slack ! 

, max := 0; row := !; col := !; 

FOR i := ! TO nosn DO 
FOR J := ! TO nodn DO 

IF (max < slackCitj]) AND (dnCj].dem> 0.00!) AND (snCiD.sup > 
THEN BEGIN 

max : = slackC i . j ] ; 
row : = i ; 
col : = J ; 

END 

ELSE IF <round(max) = round< s lackC i . j 3 ) ) AND 

<dnCj].dem > 0.00!) AND (snCi].sup > 0.00!) 

THEN IF NOT ass i gnedC i . j 3 

THEN IF <i=row) OR (j=col) 

THEN braketie 
ELSE BEGIN 

row : = i ; 
col := j; 

END; 

END; C procedure max^slack } 

/ 


0 . 001 > 



FUNCTION demand_satisf led : BOOLEAN; 

VAR 1 : INTEGER; 

BEGIN 

dQmand_8at 1 sf 1 ed := TRUE; 

FOR 1 : = 1 TO nodn DO 

IF dnC I ] . dem > 0 . 001 
THEN demand_sat 1 sf 1 ed := FALSE; 

END; 

PROCEDURE pr 1 nt_ass 1 gnment ; 

VAR i,j : INTEGER; 

BEGIN 

FOR 1 : = 1 TO nosn DO 

BEGIN 

FOR j : = 1 TO nodn DO 
WRITE (ass 1 gnment C 1 , j ]: 4 ! 0 ) : 

WRITELN; 

END; 

READLN ; 

END; 

BEGIN C Procedure F i nd_ass i gnment } 

FOR ip := 1 TO nosn DO 
FOR jp := 1 TO nodn DO 
BEGIN 

aasignmentCip.jp] := 0; 
assignedCip.jp] := FALSE; 

END; 

count := 0; 

REPEAT 

count ;= count + 1: 

Max_s lack ( ro . co) ; 

small := smal ler(snCro] .sup.dnCco] .dem > 
snCro].sup := snCrol.sup “ small: 
dnCcol.dem := dnCco].dem - small; 

Ass 1 gnment C ro . co] ;= small; 
ass i gnedCro . co] := TRUE; 

£ pr i nt__ass i gnment ; 1 

UNTIL (demand_satisf led) OR (count >= 200); 



can * * t be satisfied **’) 


IF count = 200 THEN WRITELN ( 'Demand 
tots lack : = 0 ; 

FC)F< ip 1 = 1 TO nos n DO 
FOR jp := 1 TO nodn DO 

totslack := totslack + slackC ip. jp]*assignmentC ip, jpD 
WRITELN ('Total slackness ; '. tots lack : ^ • 1 ) • 

ENO: C Procedure F i nd_ass i gnment } 


PROCEOURE Find__short_assign < nodn.nosn : nodetype; 

^ ^ • 3-ss_dem_nd; sn : ass sup nd; 

short : sl_ass_typ: 

VAR assignment : sl_ass_typ>; 

VAR ip. IP. ro.co, count : INTEGER; 

assigned : ARRAY C 1 . . nodmax , 1 . . nodmax ] OF BOOLEAN; 
totshort, small : REAL; 

PROCEDURE Min_short (VAR row, col: INTEGER); 

VAR i,j : INTEGER; min : REAL; 

PROCEDURE braketie; 

VAR penaltyl. penaity2 : REAL; 

FUNCTION next_min_in_row (iro : INTEGER) : REAL; 

VAR k : INTEGER; nextmin : REAL; 

BEGIN 

nextmin := 1000; 

FOR k := 1 TO nodn DO 

IF (nextmin > short C i ro . k] ) AND ( short C i ro , k ] > min) AND 
NOT ass i gnedC 1 ro . kl 

THEN nextmin r= shortC i ro.kl ; 
next__mi n_i n_row := nextmin; 

END; 

FUNCTION next_mi n_i n_co 1 (ico : INTEGER) : REAL; 

VAR k ; INTEGER; nextmin : REAL; 

BEGIN 

nextmin := 1000; 

FOR k :* 1 TO nosn DO 

IF (nextmin > shortCk. icol ) AND (shortCk . i co] > min) 

AND NOT ass 1 gnedCk . 1 co] 


THEN nextmin shortCk, i co] ; 
next min in col := nextmin; 


END; 

BEGIN C braketie > 
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IF j = col 
THEN BEGIN 

penaltyl := next_mi n__i n_row( i ) - min; 
penaltyZ := nGxt_min_in_row(row> - min; 

IF penaltyl > penaltyZ 
THEN row : = i ; 

END; 

I F i = row 
THEN BEGIN 

penaltyl == next_min_in_col ( j ) - min; 
penaltyZ := next_min_i n_col (col ) - min; 

If- penaltyl > penaltyZ 
T HEN col : = J ; 

END; 

END; C braketie > 

BEGIN C procedure Min_short } 

min := 1000; row := 1; col := 1; 

FOR 1 : = 1 TO nosn DO 
FOR J : * 1 TO nodn DO 

IF (min > shortCi.j]) AND (dnCj].dem> 0.001) 

AND (snCiU.sup > 0.001) C preferably check time feasibility } 
THEN BEGIN 

min : = shortC i , j 3 ; 
row : = i : 
col : = J : 

END 

ELSE IF <round(min) = round(short[ i . j 3 ) ) AND 

<dnC]3.dem > 0.001) AND (sn[i3.sup > 0.001) 
THEN IF NOT ass i gnedC i , j 3 

THEN IF (i=row) OR (j=col) 

THEN braketie 
ELSE BEGIN 

row : = i ; 
col := j; 

ENO; 

ENO; c procedure min short 3 


FUNCTION demand satisfied ; BOOLEAN: 



ow 

VAR 1 : INTEGER; 

BEGIN 

domatid_sat 1 sf i ed := TRUE; 

FOR 1 := I TO nodn 00 

IF dnCiT.dem > 0.001 
THEN demand_satisf ied := FALSE; 

END: 

PROCEDURE pr 1 nt_ass 1 gnment : 

VAR i.j : INTEGER; 

BEGIN 

FOR 1 := 1 TO nosn 00 

BEGIN 

FOR j := 1 TO nodn 00 
WRITE <ass ignmentC 1 , ) ] : 4: 0) ; 

WRITELN; 

END; 

REAOLN; 

END; 

BEGIN C Procedure F i nd_short_ass i gn > 

FOR ip 1 = 1 TO nosn 00 
FOR jp ;= 1 TO nodn DO 
BEGIN 

ass ignmentCip. jp] ;= 0; 
ass ignedC ip. Jp] != FALSE; 

END; 

count j= 0; 

REPEAT 

count ; = count + 1 : 

Mi n_short (ro , co) ; 

small := smal lerlsnCro] .sup.dnCco] .dem >; 

snCro].sup snCroD.sup ~ small; 
dnCcoJ.dem := dnCcoj.dem - small; 

Ass ignmantCro.co] := small; 
ass i gnedCro , co] := TRUE; 

C pri nt_ass ignment ; > 

UNTIL <demand_satisf led) OR (count >= 200); 

IF count = 200 THEN WRITELN < 'Demand can"t be satisfied **' > 


totshort ;» 0; 
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FOR IP : = 1 TO nosn DO 
FOR jp := 1 TO nodn DO 

totshort .= totshort f shor t [ i p . ,p3*ass . gnment C i p . j p] ; 
WRITELN ('Total shortness : totshort : ^ : 1 ) : 

END; C F i nd_8hort_as3 1 gn } 

PROCEDURE f ind_short_ex_aray ( s. t : nodetype; ex_arc : arcset; 

VAR short_aray : arcarr; VAR stackptr : INTEGER ); 

LABEL 20; 

VAR maxdist, temp_dist : REAL; 

Nset , temp, perm : nodeset; 
shordist : ARRAY Cl..nodmax] OF REAL; 
i , j , arc_no : INTEGER; 
store, storl : nodetype; 

FUNCTION min_temp_di st_nod ; nodetype ; 

VAR minimum : REAL; k ; nodetype; 

BEGIN 

minimum := maxdist; 

FOR k i« 1 TO no_nod DO 
IF NOT <k IN perm) 

THEN IF shordistCk] < minimum 
THEN BEGIN 

minimum ;= shordistCk]; 
mi n_temp_d i st_nod ;= k; 

END; 

END; C min;__t 0 mp_dist_nod ) 

BEGIN C Procedure find short aray ) 

Nset 1 = C 1 . • no__nod3 ; 

perm : “ Cs]; temp Nset - Cs3; 

IF NOT <<s IN Nset) AND (t IN Nset)) 

THEN BEGIN 

WRITELN ('Node out of network data range !!'); 

GOTO JOO; 

END; 

C Set shortest d.st for other nodes = sum of all dist.s in network 
maxdist 0; 


FOR i ; = 1 TO 2 *no_arc DO 



maxdist := maxdist + arc[i].dist: 

FOR 1 : = 1 TO no__nod 00 

shordistCi] •= maxdistj 
shordistCs] := 0; 
o tore : = s ; 

WHILE NOT (t IN perm) 00 
BEGIN 

FOR j:= pointCstore] TO po i nt [ store+1] -1 00 
BEGIN C Get the node in the ad)uscency list of node ’store' 
called storl ' , update its shordist if required } 
arc_no := ftraceCj]; 

IF not<arc_no IN ex_arc> 

THEN temp_dist := shordistCstore] + arcCarc_no3 . d i s 
ELSE temp_dist := 10000; 
storl := arcCarc_no] . nodeZ : 

IF storl IN temp 

THEN shordisttstorl ] := smaller (shordistCstorl I , 

t emp_d 1 s t > ; 

END; C For adjuscent nodes update shordist } 
store := min_temp_dist_nod; 
temp := temp - [store]; 

• perm ;= perm + [store]; 

END; C while not (t in perm ) DO } 

C find the arcs in the shortest path } 

perm := [t]; temp := Nset - [t]; 

* 

store := t; stackptr ;= 0; 

I WHILE NOT (8 IN perm) 00 
BEGIN 

FOR point[8tore] TO point[store+l ] -1 DO 

BEGIN 

arc_no := rtrace[)]; 

I IF <shordist[ arc[arc_no] . nodel ] = shordi st[store] 

- arc [arc_no] .diet) AND NOT (arc_no IN ex_arc) 

THEN BEGIN 

store := arc[arc_no] .nodel ; 

GOTO 20; 

END; 

END; 
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stackptr ; ■ stackptr + 1: 
short aray[8tackptr] ■• = arc_no; 



temp ;= temp - [store]; 
perm := perm + [store]; 

END; C while not <s in perm ) do ) 

END: C Procedure find short ex aray J 

procedure Update_deiay_info<dele : REAL; flo : REAL; 

£ 

VAR del_aray : del_type; VAR wt_dele : wtd delay typ; 

VAR stackptr : INTEGER); 

VAR int_dal : INTEGER; 

BEGIN 

int_del := round<dele); 

IF int_dei > stackptr 
THEN stackptr ;= int_dei: 

de l_aray C 1 nt_de 1 ] := del_arayC int_del ] + round<fIo>; 
wt_deleC int_del ] := wt_deleC int_del ] + dele*flo: 

END; 


BEGIN C mam > 

WRITE ('Do you want to begin with infinite capacity of arcs? ') 
READLN Cans); 

IF (ans = 'y ' ) OR (ans = 'Y' ) 

THEN RESET ( i np . ' i npsup . dat ' > 

ELSE RESET < inp, ' inpsupcap.d’ > ; 

READLN (inp,no_nod. no_arc); 

FOR 1 i I = 1 TO no_arc 00 

WITH arcC i i 3 DO 

IF (ans =» 'y'> OR (ans = 'Y') 

THEN READLN (inp, nodel . nodeZ.dist) 

ELSE READLN (inp. nodel. node2.dist.capa); 

FOR 11 : = 1 TO no_arc DO 

WITH arcC i i 3 DO 

BEGIN 

cost_9 ;= dist / 10; 
cost w cost_g /1. 5; 
dist := dlst / (?0*4); 

IF (ans = 'y'> OR (ans = 'Y') 


/ 


THEN capa ;= 1000; 
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ENOt 

FOR II : = no_arc + l TO 2*no_arc DO 
BEGIN 

arcC 1 1 ] . nodel := arc[ i i -no_arc] . nodeZ ; 
arcC 1 1 ] . nodeZ := arcC i i -no_arc] . nodel ; 
ai'cC i 1 3 . cost_q := arc[ i i -no_arc 3 . cost_g ; 
arcC 1 1 3 . cost_w := arcC i i -no_arc3 . cost_w; 
arcCii3.dist := arc[ i i -no_arc3 . dist ; 
arcCii3.capa := arcC i i -no_arc3 . capa ; 

END; 

WRITE ('Give period length : ' ) ; READLN (period_lenth> ; 
WRITE (’Give time horrizon : *>; READLN (T_hor); 

IF T_hor < period_lenth 
THEN BEGIN 

WRITELN ('Time horrizon must be > period lenth ’) 
GOTO Z5; 

END; 

FOR 11 : = 1 TO 2*no_arc DO 

WITH arcCii3 DO C initialise 1 
FOR tt := 1 TO T_hormax DO 
BEGIN 

res_capCtt3 ;= capa; 
loadCtt3 := 0; 

END; 

POR II : = 1 TO no_nod DO 
WITH nodeCii] DO C Initialise 3 
FOR tt := 1 TO T_hor DO 
BEGIN 

dema_supl Ctt3 Ot 

requ_avai [tt3 := 0: 

END; 

WHILE NOT EOF (inp) DO 
BEGIN 

READ <ihp. nod.ds.time) : 

WITH nodeCnod3 DO 


dema_9UplCtime*63 = = ds 

IF EOLN ( inp) 

THEN READLNdnp)! 


END; 
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RESET < J np , ' 1 npdem. dat ' ) ; 

WHILE NOT EOF < mp) 00 
BEGIN 

READ (inp, nod. ds, time) j 
WITH nodeCnod] 00 
dema_supl C t ime*63 := -ds; 

IF EOLN <inp) 

THEN READLN(inp); 

END; 

RESET (inp. • inprowag.daf); £ Check and delet to read i 

WHILE NOT EOF (inp) 00 

BEGIN 

READLN C 1 np , nod , ds , t ime) ; 

WITH nodeCnodT 00 

requ_avai Ctime*6] := requ_avai C t ime*6 ] + ds; 

END; 

CLOSE (inp); 

FOR i 1 : = 1 TO no_nod 00 

WITH nodeC i i ] DO 

cum_dem_sup_g i= 0; 

FOR 1 i : = 1 TO no_nod 00 
degreeC i i ] ; = 0 ; 

FOR i i ! = 1 TO 2*no_arc DO 
WITH arcC i i 3 DO 

dQgreeCnodeZ3 := degree[node23 + 1; 
cumsum ; = 1 ; 

FOR 11 ! «= I TO no_nod+-l DO 

BEGIN 

pointCii] t= cumsum; 

cumsum ;= cumsum + degreeCii3: 

END; 

FOR 11 ;= 1 TO no_nod+l DO C Ini t lal i sat i on > 

BEGIN 

fmovptrCii3 := point[ii3; 
rmovptrCii3 s= pointCii3; 


END; 

FOR 11 ; s 1 TO 2*no_arc DO 

WITH arcC 1 1 3 00 

BEGIN 


t once } 


/ 



ftraceC fmovptrCnodel ] ] .= 
fmovptrCnodel] := fmovptr [nodel ] + i, 
rtraceC rmovptrCnode2] ] ,= ij. 
rmovptrCnodaZ] s= rmovptrCnode 2 ] + l, 

END? 

FOR 11 : = 0 TO T_hor DO 

BEGIN 

de 1 ay_f reqC 1 1 ] := q. 
wtd^de lay [ I i ] •= O; 

END; 

REWRITE <out2» * outroute . dat * > I 
REWRITE < ou t 5 , * outcost . dat ’ ) ; 

WRITELN <out^, 'Period Goods Cost Empty Wagon Cost’); 
par 1 od_no i =« 1 ; 

WRITELN; WRITELN ('Goods wagon scheduling'); 

WRITELN (’Period no : ' ,period_no: 2) ; 

T_ t= period_no* period_lanth: 

FOR 1 i s = 1 TO no_nod DO 
WITH nodeC i i ] DO 
BEGIN 

FOR tt :» (pariod_no-l )*period_lenth + 1 

to period_no*per iod_Ienth DO 
cum_dem_sup_g := cum_dem_sup_g + dema_suplCtt3 ; 
IF cum_dem_8up_g < -0.001 
THEN dem_nod_8et ;= dem_nod_set + Ciil 

END; 

mm no_,nod + no_arc * 2; C No. of constraints ) 
ni {“ 2»no__arc; C No of decision variables > 

FOR 11 ; “ 1 TO 2*no__arc DO 

cmCrnm+iil :* arcC i i ] . cost_g; 

FOR II I » 1 TO mm DO 
FOR jj ;= mm+1 TO mm+ni DO 
AmC i 1 . j j 3 t = 0 ; 

FOR i i ; = 1 TO no_nod DO 
WITH nodeCii] DO 
BEGIN 

£ For all the arcs incoming at node ii > 

FOR jj pointCiil TO pointCii+13 -1 DO 


BEGIN 


87 


kk : = mm + rtraceCj)]; 

IF cum_dem_sup_ g <= o 
f HEN AmC J 1 , kk] : = 1 
ELSE AmCii.kk] := -!• 

END! 

C For all the outgoing arcs from node ii } 
FOR ij := pointCii] TO pointCii+1] -1 DO 
BEGIN 

kk := mm + ftraceCjj]; 

IF cum_dem_sup_g <= 0 
THEN AmC i i .kk] -1 
ELSE AmC i i . kk] : = 1 ; 


END: 

IF cum_dem_sup_g < -0.001 
THEN 1 nequal 1 tymC i 1 ] := ’>’ 

ELSE IF abs ( cum_dem_sup_g) <= 0.001 
THEN 1 nequal 1 tymC 1 1 ] := ’=' 
ELSE inequal 1 tymC 1 1 ] := ' < ' ; 
bmiCii] := abs<nodeC 1 1 ] .cum_dem_sup_g) ; 
END; C for 11 : = 1 TO no_nod do 1 

WRITELN (‘Node cum_dem_sup_9 ' > : 

FOR 11 ! = 1 TO no_nod DO 

WRITELN ( i i ; 4, nodeC i i ] . cum_dem_sup_g : 8 : 1 > ; 

Ji := 0; 

FOR 11 := 1 TO no_arc*2 DO 

BEGIN 


j j } j + 1 5 

AmC i 1 +no_nod , rw)+ j j ] := 1: 

1 nequal i tymC i 1 +no__nod] ;= 
bmi Cno_nod+ 1 1 3 arcC 1 1 3 • capa; 


END; 

simplex (mrn.ni .Cm, Am. inequal itym.bmi 
pariod_cost s= 0; 

IF mincoat__feasible THEN 

FOR 11 -•= 1 TO 2*no_arc 00 

period_cost t= period_cost + 
snn ;= 0; C No.of sup nodes > 
dnn ;= 0; t No.of dem nodes 1 


bmo,arc_f lo_aet.mincost_fea3ible) 


bmoC i i 3*arcC i i 3 - cost_g: 


FOR ii : = 1 TO no nod+l DO 


RPr tM 



ddegL i i J : = 0 ; 
sdegC 1 1 ] : = 0 ; 


END; 

FOR i 1 : = I TO no_nod DO 

WITH nodeC i i ] DO 
BEGIN 

FOR tt := 1 TO T_ DO 
IF dema_supl Ctt] < -0.001 
THEN BEGIN 

dnn : = dnn + 1 ; 

ddegC 1 1 ] := ddegCii] + Ij 

WITH dmnCdnn] DO 

BEGIN nd := n; 

dam != abs (dema_supl C t t 3 ) j 
tim := tt; 

END; 

END 

ELSE IF dema_suplCtt3 > 0.001 
THEN BEGIN 

snn : = snn + 1 ; 
sdegCii] := sdeg[ii3 + 1; 

WITH spnCsnn] DO 
BEGIN nd := i i ; 

sup != dema_supl Ctt] ; 
t im : = tt : 

END: 

END: 

END; 

IF mi ncost_f eas 1 ble 
THEN BEGIN 

FOR li != 1 TO snn DO 
FOR J J : = 1 TO dnn DO 
BEGIN 

i no spnC i 1 3 • ; 

J no := dmnCjj3.nd; 

slakCii.Jj3 := dmnCnl.t.m - spnCiiD-tm 

- shortest_dist< ino, jno) 


END: 

f i nd_as s i gnment ( dnn . snn . dmn . spn . s 1 ak , ass i gn ) 


END 
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ELSE BEGIN 

FOR it ! = 1 TO snn DO 
FOR j j r = 1 TO dnn DO 
BEGIN 

ino != spnC i i ] .nd; 
i no ! = dmnt ) j ] . nd ; 

slakCii.jj] := shor tes t_d I 3 1 ( 1 no , j no) t 

END; 

F t nd_short_ass i gn{dnn , snn , dmn , spn , s lak , ass i gn> ; 
FOR I i := 1 TO 2*no_arc DO 
WITH arcC i i ] DO 

FOR tt := ( pe r t od_no- 1 ) #per i od_l enth+1 TO 
per I od_no»per 1 od_l enth DO 
res_cap[tt] := capa/6i 

END; 

FOR II := 1 TO 2*no_arc DO 

WI TFi arcC i i 3 DO 

BEGIN 

tmp_capa ; = bmoC it]; 

IF tfnp_capa < 0.001 
THEN deleted := TRUE 
EL.SE deleted := FALSE; 

END I 


cumsum ; = 1 ; 

FOR ii ; = 1 TO no_nod+l DO 
BEGIN 

dptCi i] s= cumsum; 

cumsum s= cumsum + ddegCii]; 

END; 

cumsum : = 1 ; 

FOR 1 i ;= 1 TO no_nod +1 DO 
BEGIN 

sptC i I 3 != cumsum; 

cumsum j= cumsum + sdcg[ii3; 

END; 

FOR i i : = 1 TO dnn DO 

WRITE (dmnC i i 3 . nd: J > ! 



WRITELN: 


WRITE <’ ->. 

FOR n : ^ 1 ro dnn DO 

WRITE (dmnC i 1 3 .dem: 5 :0> ; 

WRITELN: 

WRITE <• •); 

FOR 11 : = 1 TO dnn DO 

WRITE <dmnC 1 1 3 . t im: > ) : 

WRITELN; 

FOR s_nod : = 1 TO no_nod DO 

FOR jj>= sptCs_nod3 TO sptCs_nod+l 3-1 DO 

BEGIN 

WRITE < spnC j j 3 . nd : 5 , spn[ j j 3 . sup : J 0 , spnC j j 3 
FOR cur_nod := 1 TO no_nod DO 
FOR kk : = dptCcur_nod3 TO dptCcur_nod+l 3-1 DO 
WfllTE (ass 1 gnC J ) , kk3 : 5 ! 0 ) ; 

WRITELN: 

END; 

WRITELN <out2, 'Goods Period no : ' , per i od_no : 2 ) ; 

WRITELN (out2.’Flow Node Time 

cn : ® 1 : s_nod : = 1 j 
IF mi nco8t_f eas ibl e THEN 

WHILE <dem_nod_set <> C3) AND (cn<400> DO 
BEGIN 

WHILE <nodeCs_nod3 . cum_dem_sup_g < 0.001) 

AND (s_nod <=* no_nod) DO 

s_nod := s_nod + 1: 
cur__nod s= s_nod; 

WRITELN Ccn ; '.cn;2,’ s_nod'. s_nod:4,' 

' , nodeC8_nod3 . cum_dem__sup_g 1 4 : 0) ; 

min_arc_capa := 1000; 

no : = 0 ; 

travl^tim s =' 0; fio_sent 0; 

FOR ii ;= 1 TO 2*no_arc DO 

arcCiil.Goter TRUE ; 

7 ?: REPEAT 

WRITELN <'cur_nod : ' . cur_nod : 4) : 


t imi > 


> ! 


sup 


WRITE ('ft 

FOR ii pointCcur_nod3 TO pointCcur_nod+l 3-1 DO 



WITH arcC f traceC i i ] 3 DO 
BEGIN 
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ft^in_arc_capa> ; 


200 ? 

t im-spnC j store 3 

inf eas ibi i i ty 
- aval l_t im; 
delay : 5 : 2 ) ; 


ft := ftraceCii]; WRITE (fti4); 

IF <ft IN arc_fio_set) AND NOT deleted AND enter 
THEN BEGIN 

WRITELN ( • enters’ ) ; 

arc_stor := ft; 

no ? = no t 1 ; 

arcs_aray[no] := ft; 

travi_tim := travl_tim + dist; 

prev_mi n_a redcap := min_arc_capa; 

min_arc_capa : = smai i er(arc[ f 1 3 . tmp_capa , 

IF node2 IN dem_nod_set 
THEN BEGIN 

ass := 0; 

FOR JJ:= sptCs_nod] TO sptCs_nod+l 3-1 00 
FOR kk!= dpt[node2] TO dptCnodeZ+1 3-1 DO 
IF ass ignC j ) .kk3 > 0.001 
THEN BEGIN 

ass := assignC j j .kk3 ; 
j store : = ) J i 
kstore ;— kk ; 

GOTO ZOOiCexit for loop! 

END; 

WRITELN ('Demand nod * . nodeZ : 4> ; 

IF ass > 0.001 
THEN BEGIN 

avail_tim dmnCkstorel. 

IF avail_tim < travl_tim 
THEN BEGIN 

WRITE ('time 

delay := travl_tim 

WRITELN < 'delay ' , 

s tim := spnCjstore 


f i m * 
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round< t rav i_t im) ; 


d tim := s tim + 


> 0-001 


END 

ELSE BEGIN 
d_tim := dmnCkstore] - t im; 

delay := 0; 
“ rounddravl tim) 


-round<travi tim) 


THEN s tim d tim 


deiay_f req, wtd^^delay ,ptr) ; 
- flo sent; 


ELSE s_tim := 0; 

END; 

flo_sent := smal i er(min_arc_capa, ass): 
Update_de lay_i nf o (de lay , f lo_sent , 

assignC j store , ks tore] := assignC jstore,kstore3 

WITH dmnCkstore] DO 
BEGIN 

dem := dem - fio^sent; 
FOR jj := i TO snn DO 
assignC|),kstore3 


assignC n .t<store3 - flo^sent; 
nodeCnd] . dema^supi C t im] f lo_sent ; 

nodeCnd] . requ^avai Cd_t im] + flo_sent; 


nodeCnd] . dema_supl Ct im] s = 
nodeCnd] . requ_avai Cd_t im] : = 
WRITELN ( n im d ’ , 


d_t im: 2 I 

END; 

WITH spnCjstore] DO 
BEGIN 

sup := sup - flo^^seni 
FOR kk := 1 TO dnn DC 
assignCjstoretkk] 


ass ignC jstore,kk] - flo_sent; 


WRITELN <’tim s *,s_tim:2>! 

nodeCnd] -dema_supl Ct im] 


nodeCnd] - dema^supl Ct im] ^ • 



THEN nodeCnd] . requ avaiCs timD := 

odeCnd’J . roqu_av.x i L !i__t im3 - flo_sQnt: ~ ~ 

END; 

CWRITE arcs aray ok fwdJ 
kk := d_tim: 
rk := kk; 

FOR j j ! = no DOWNTO 1 DO 
WITH arcC arcs_aray [ j ) ] 3DO 
BEG-IN 

FOR tt := kk DOWNTO 

kk-round(dj st> DO 


res__capC 1 1 3 - f 1 o_sent ; 
loadCtt3 + fio sent; 


arcs_arayC j j 3 : 2 > ; 


re8__capC 1 1 3 ; 5 I 0 . ' Load 


BEGIN 

r8S_capCtt3 := 


loadCtt3 := 

I 


* , ioadCtt]i?:0)? 


END; 


END; 

WRITELN (*Arc no 

1 

I 

WRITELN <*Res capa 

WRITE ( 'rk : * ,rk: 3:0> :1 
rk != rk - dist; 1 

WRITELN(rk:3:0) ; 
kk : = round! rk) ; 


WRITELNCflo sent ’ . f lo_sent : 4 : 0 > ; 
END C ass > 0.001 } 

ELSE GOTO 250; £ ass <= 0.001 } 


END 

ELSE GOTO 250;Cnode2 not in dem_nod_set} 
nodeCs_nod3 . cum_dem_sup_g := 

nodeC 8__nod3 , cum__dem_sup_g "• flo_sent; 

FOR JJ := 1 TO no DO 
WITH arcCarcs_arayC ) J 33 DO 
BEGIN 

tmp_capa := tmp_capa - flo_sent; 
IF tmp_capa < 0.001 
THEN deleted := TRUE; 
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nodeCnodeZ 3 . cuni_defn_sup_g : = 

nodeCnode?. J - cum_dem_8up_g + flo_sent! 

IF nod0[node2] .cum_,dem .sup_q > -0.001 
I HEN dem_nod_set := dem_nod_set - CnodeZ]; 
WRITELN ( ‘flag 2 * ) ; 

GOTO 100: 

END; CIF ft in arc_flo_st & not arcC f t ] . de le tdJ 
E.NO: Cfor ii := fpoint.. & WITH arc[ii3 do> 

WRITELN; 

{s_nod had surplus suplies which could not be absorbed 
WITH current assignments) 

IF s_nod <= no_nod 
THEN IF no > 0 
THEN BEGIN 

arctarcs_arayCnoD 3 . enter := FALSE; 
min_arc_capa := prev_mi n_arc_cap ; 
cur_nod := arcCarcs_aray[no3 3 . nodel ; 
no : = no - 1 : 

WRITELN ( 'Backtracking' ) : 

GOTO 75: 

END 

ELSE BEGIN 

9 nod : - 9_nod + 1 ; 

WRITELN ('Futile search. 

s^nod Quit., surplus assumed') 

END 

ELSE BEGIN C s_nod > no_nod ) 

WRITELN ('Failure with current assignments, 

discarding them' ) i 

8__nod ; * 1 ; 

FOR 11 :® 1 

FOR jj := 1 TO dnn DO 

IF (spnCii3.sup>0.001) AND (dmnC JJ 3 • dem>0 . 001 ) 
THEN ass ignCi 1 .) J3 /•= smaller 

(spnC 1 1 3 . sup.dmnC ) } 3 -dem) ; 

END: 

cn :=» cn + 1? 
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100 : 


5 ? 0 : 


575 : 


cn r = cn 4- 1 ; 

UNTIL (flo_sent > 0.001) OR (cn>400>; 

WRITELN < ' f iag5 • ) : 
rk ! = s_t im? 

IF flo_sent > 0.001 
THEN BEGIN 

WRITE (out2 , f lo_sent :6 : 2 , ' *>5 

FOR 1 1 : = 1 TO no DO 

WITH arcCarcs_arayC i i 1 3 DO 

BEGIN 

WR ITE (out2,nodel:8,rk:8i2>i 
rk f = rk + d i s t ; 

END? 

IF no <> 0 

THEN WITH arc [ arcs_aray C no 3 3 DO 

WRITELN (out2,node2:8,rk:8i2>; 

END? 

END C while dem_nod_set <> C3 3 

ELSE C mincost infeasible > 

WHILE (dem_nod_set <> [3) AND (cn<180) DO 
BEGIN 

WRITE (*Dem nod set : *)? 

FOR 11 : = 1 TO no_nod DO 

IF ii IN dem_nod_set 
THEN WRITE (ii:2); WRITELN; 

FOR 1 i ; = 1 TO snn DO 
FOR j ) ! = 1 TO dnn DO 

IF (ass 1 gnC i i , j j 3 > 0) and (dmn[jj3.nd IN dem_nod__set ) 
THEN GOTO 550; 
s_nod j= spnCii3.nd; 
d__nod := dmnCjj3.nd; 

avail_tim := dmnCjj3.tim - spnCii3.tim; 

flo_sent != smaller < spnC i i 3 . sup. dmn[ j j 3 . dem) ; 

flo_sent ;= smaller (flo_sent. ass i gnC i i * j j 3 > ; 

istore := ii; jstore := jj; 

rk ! = dmnC | j 3 . t i m? 

exci_arc i= C3i tried := FALSE? 

f i nd_short_ex_aray ( s_nod, d_nod . exc I_arc .arcs_aray , no) t 
C Here arcs_aray*s storage is in reverse order > 
mot ion t i m : = 0 ; 
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580f 


FOR i 1 f = 1 TO no DO 

motion_tim != motion_tim T arc[arc9_aray C i i 3 3 . d i b t j 
IF (nnotion_tim > avaii_tim) AND NOT tried 
THEN BEGIN 

WRITELN <*Time i nf eas i bl e * > ; 
tried s= TRUE; 

IF exci_arc <> C3 
THEN BEGIN 

e X c 1 _a r c : = C 3 ; 
goto 
END; 

END; 

WR I TELN < * s_nod ' , s_nod : 5 , ' d_nod * . d_nod : 5 > ; 

WRITE < 'Short aray*); 

FOR i 1 ! = no DOWNTO 1 DO 

WRITE (arcs_aray [ i i 3 : 4) : 

WRITELN; 

C We check the residual capacities of the arcs, and time feasi} 
travl_tim ;= 0; 

FOR i 1 := 1 TO no DO 

WITH arc[arcs_aray C i i 3 3 DO 

BEGIN 

tim_dt 1 = rk~ t rav l_t imr 
t i m_s t s = t i m_d t - diet; 

IF tim_st < 0 THEN tim_st ;= 0; 
passed ;= TRUE; 

FOR tt round! t im_dt ) DOWNTO round! t i m_8 t ) DO 

BEGIN 

IF res_capCtt3 <= 0 
THEN BEGIN 

tim_dt tt - 1; 
t i m_s t : = t i m_d t d i s t ; 
passed := FALSE; 

END 

ELSE flo^sent != sma 1 1 e r ! f 1 o^sen t , res^capC t t 3 ) i 

END; 

IF NOT passed 

THEN IF degree[arc[arcs_aray C i i 3 3 . node23 > 1 

THEN excl_arc := excl_arc + [arcs_arayC I i 3 3 
ELSE tried ;= TRUE; 
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tra victim j= rk ~ tim^sts 

IF NOT passed AND ( < t rav l__t im<ava i l_t im) OR tried) 
THEN GOTO >80 i 

slack^tim s= avail_tim - travl_timj 
IF (slack_tim < 0) AND (mot j on_t i m<ava i l_t im) 

AND NOT tried 

THEN BEGIN 

WRITE ( ’Excl Arc ’ ) i 
FOR jj 5= 1 TO 2*no_arc DO 
IF j j IN excl_arc 
THEN WRITE ( j J : 4) : 

WRITELN; 

WRITELN (’Time infeasible (temp)’); 

GOTO >7>; 

END; 

C else if motion_tim > tim then delays > 

END; 

IF avail^tim < travl_tim 
THEN BEGIN 

WRITE (’Time i n f eas i b i 1 i ty ’); 
delay ;= travl_tim “ avail_tim; 

END 

ELSE delay 0; 

WRITELN (’delay ’ , de 1 ay r > : 2 ) ; 

Updat e_de 1 ay_i nf o ( de 1 ay , f 1 o_sent , de lay_f req , wtd_de lay , pt r ) t 
FOR i i ' ; = no DOWNTO 1 DO 
WITH arcCarcs_aray [ i 1 ] 3 DO 
BEGIN 

FOR tt := round( t im_s t ) TO round! t im_dt ) DO 
BEGIN 

re3__capCtt] := res_cap[tt3 "* flo^sent; 
load[tt3 1 = loadCtt] -F flo__sent; 

END ; 

WRITELN (’arc no ’ , arcs_aray C i i 3 * 2 , 

’ Res capa ’ , re8__capC 1 1 3 : 5 s 0 , ’ load ’ , 1 oadC 1 1 3 : > : 0 > ; 

WRITE Ctims t t ’ , t i m_s t ; > : 0 ) ; 

WRITELN (tim_dt :>:0) ; 

’ period_cost ?= period_cost + flo^sent* cost^g; 


END 
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WRITELN (*fio sent * , f i o_sent r 4 : 1 ) ? 

assignCistore,|store3j= assignCistore.jstore] "fl o^sent i 
spnC i s tore 3 . sup := spnC i s tore 3 . sup ~ fio_sent; 
time := spnC i store 3 - t.im; 

WITH nodeCs_nod3 DO 
BEGIN 

dema_sup 1 C t I me 3 := dema_sup 1 C t i me 3 “ flo_sent| 
requ_ava i C t ime 3 := requ_ava i C 1 1 me 3 ** flo_senti 
cum__dem_sup_g := cum_dem_sup_g - flo__sent; 

END? 

dmnC j s tore 3 - dem := dmnC j s tore 3 - dem - flo__sent; 
time != dmnC j store 3 . t im; 

WITH nodeCd_nod3 DO 
BEGIN 

dema_3up 1 C t ime 3 dema_3upl C t ime 3 + flo^sentj 
requ_ava i C t ime 3 := requ_ava i C t ime3 ^ flo^sentj 
cum_dem_sup_g := cum__dem_sup_g -t flo^sentj 
IF cum_dem_sup_g > -0.001 

THEN dem_nod_set := dem_nod_set - Cd_nod3: 

END? 

IF flo_sent > 0.001 
THEN BEGIN 

WRITE <out2. f lo_sent:6i2, ’ •>? 

FOR i i : = no DOWNTO 1 DO 
WITH arcCarc9_aray C i I 3 3 DO 

WRITE CoutZ.nodel :8. tim_st :8i£>i 
IF no <> 0 

THEN WITH arcCarc3_aray C 1 3 3 DO 

WRITELN(out£.node2 : 8, ( t im_s t-td i st ) : 8 s £ > 

ELSE WRITELN ( out 2 . 9 _nod : 8 , spnC i store 3 . t imx 6 . d_nod 1 8 . 
dmnC jstore3.tim:6)? 


END? 

cn J = cn tl ? 


END? C while <dem_nod_set <> C3> and (cn<180) do) 

IF dem_nod_set <> C3 THEN 

WRITELN (’Exiting while loop !! goods demand unsatisfied’)? 
WRITE (out^,period_no:4, period_cost : 12?2> ? 
total cost 1 = totai_cost + period^cost? 


300 ? 
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C Empty wagon scheduling 

y 

WRITELN; WRITELN < 'Empty wagon scheduling’); 
WRITELN (’Period no : ’ , per i od_no : 2 > ; 

period_cost := 0; 
dem_nod_set := C3; 

FOR 11 := i TO no_nod DO 

WITH nodeC 1 i 3 00 

BEGIN 

IF period_no = 1 
THEN BEGIN 


cum_dem_sup_w := 0; 

FOR tt := 1 TO period_lenth DO 
IF requ_ava i C t t 3 > 0.001 

THEN C:um_dem_sup_w ;= cum_dem_sup_w + requ_ava i C tt 3 

END 

ELSE C period_no > 1 } 

FOR tt := ( per 1 od_no~l ) *per i od_l enth + 1 
to per 1 od_no*per 1 od_l enth DO 
C Consider demand for wagons in current time period > 

BEGIN 

IF requ_ava i C t t 3 < -0.001 
THEN BEGIN 


+ requ^^avai C tt 3 ; 


cum_dem_sup_w : = cum_dem_sup__w 

dem_nod_set : = dem_nod_set + C i i 3 ; 
END; 


C Consider suply in wagons in previous time period > 
IF requ_avai C t t 3 > 0.001 
THEN cum_dem_sup_w := cum_dem_sup_w 


+ requ_avai Ctt 3 ? 


END; 


END; 

IF period_no = 1 
THEN BEGIN 

WRITELN (’Empty wagon scheduling is assumed 
to be done in the previous time horizon*); 

WRITELN (out>) ; 


C Also mention the wagons for which assume this > 
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GOTO 600; 

END; 

WRITELN ('Node cum_dem_sup__w* ) ; 

FOR i i : = 1 TO no_nod 00 

WRITELN < i i : 4 , nodeC i i 3 . cum_dem__sup_w : 8:1); 

WRITE ( *Dem nod set : ' ) ; 

FOR 11 : = 1 TO no_nod DO 

IF 11 IN dem_nod__set 
THEN WRITE <ii:2>; WRITELN; 
dnn := 0; C No. of dem nodes > 
snn := 0; C No.of sup nodes > 

FOR 11 : = 1 TO no_nod+l DO 

BEGIN 

ddegC i i 3 : = 0 ; 

sdegC i i 3 : = 0 ; 

END; 

FOR 11 ; = 1 TO no_nod DO 

WITH nodeCi i3 DO 
BEGIN 

FOR tt != per i od_i enthfl TO T_ DO 
IF requ_avai C tt 3 < -0.001 
THEN BEGIN 

dnn : = dnn + 1 ; 

ddegC i 1 3 := ddegC i i 3 + 1; 

WITH dmnCdnnl DO 
BEGIN nd := ii; 

dem := abs< requ_ava 1 C t t 3 > ; 
t im : = t t ; 

END; 

END 

END; 

FOR i i : = 1 TO no_nod DO 

WITH nodeCi i3 DO 

BEGIN 

FOR tt := 1 TO T_ DO 
IF requ_avai [tt 3 > 0.001 
THEN BEGIN 


snn : = snn + 1 ; 
sdegCiil sdegC i i 3 + 1; 
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END; 
FOR i 1 
FOR j 3 
BEGIN 


WITH spnCsnn] DO 
BEGIN nd := ii; 

sup := requ_ava 1 C tt 3 ; 

t I nri : = 1 1 ; 


END; 

END; 


1 TO snn DO 
1 TO dnn DO 


1 no : = spnC 1 1 3 . nd ; 
j no : = dmnC j j 3 . nd ; 

slakCii.jj] := shortest_di st ( i no , jno) ; 

END; 


f i nd_shor t_ass i gn < dnn , snn , dmn . spn , s 1 ak , ass i gn > ; 
C WRITE dem/sup nod set ok 3 
cumsum : = 1 ; 

FOR ii := 1 TO no_nod+l DO 
BEGIN 

dptCii3 := cumsum; 

cumsum := cumsum + ddegC i i 3 ; 

END; 

cumsum : = 1 ; 

FOR 11 ; = 1 TO no_nod +1 DO 

BEGIN 

sptCiil r= cumsum; 

cumsum ;= cumsum + sdegCiil; 

END; 

WRITELN (Mi dpt spt ’ > ; 

FOR 11 ;= 1 TO no_nod+l DO 

WRITELN (ii:3, dptCii3:3, sptCii3:3); 

WRITELN ( 'Assignments ’ > ; 

WRITE (' 

FOR 11 := 1 TO dnn DO 

WR I TE < dmn C i i 3 - nd : 5 > ; 

WRITELN; 

WRITE ( ’ M ; 


FOR 11 := 1 TO dnn DO 



5 ? 0 : 


> 7 >: 
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WRITE (dmnC i 1 3 .dem: 5 :0> ; 

WRITELN; 

WRITE <* *>; 

FOR 11 : = 1 TO dnn DO 

WRITE (dmnCi i].tim:5>; 

WRITELN; 

FOR s_nod : = 1 TO no_nod DO 

FOR jj:= sptCs_nod] TO spt C s_nod+l 3 - 1 DO 

BEGIN 

WRITE <spnCj]3.nd:5, spnCjj3.sup:5:0,spnCjj3.tim:5); 

FOR cur_nod : = 1 TO no_nod DO 

FOR kk : = dpt Ccur__nod3 TO dpt C cur_nod+l 3 “ 1 DO 

WRITE <ass i gnC j j , kk3 : 5 : 0 > ; 

CWR I TE ( short est_d i s t ( spnC ) j 3 . nd , dmnCkk 3 . nd) : 5 : 0 > ; > 

WRITELN; 

END; 

WRITELN (out2 .* Empty Wagon Movement Period no : * , per i od_no : 2 ) 

WRITELN (out2,*Flow Node Time *>; 

cn : = 1 ; s_nod : - 1 ; 

WHILE <dem_nod_set <> C3> AND (cn<200> DO 
BEGIN 

WRITE (’Dem nod set : ’>; 

FOR 11 ; = 1 TO no_nod DO 

IF 11 IN dem_nod_set 
THEN WRITE <ii:2>; WRITELN; 

FOR 11 : = 1 TO snn DO 
FOR j j : = 1 TO dnn DO 

IF ass 1 gnC 1 1 , j } 3 > 0.001 
THEN GOTO ?50; 
s_nod ;= spnCii3.nd; 
d_nod ;= dmnCjj3.nd; 

avail^^tim := dmnC)|].tim - spnCii3*tim; 
flo_sent ;= smaller ( spnC i i 3 . sup , dmnC j j 3 . dem> ; 
flo^sent ;= smaller <flo_sent, ass i gnC i i , j j 3 ) ; 
istore := ii; jstore t= jj; 
rk : = spnC i i 3 . t im; 
excl_arc C3; tried FALSE; 

find short_ex_aray <s_nod,d_nod, excl_arc ,arcs_aray ,no> ; 


C Here arcs_aray’s storage is in reverse order > 
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motion_tim := Oj 
FOR 11 := 1 TO no DO 

:= motion_tim + arcCarcs_aray C i i 3 ] . di st ; 
IF <mot ion__t im > avaii_tim> AND NOT tried 
THEN BEGIN 

WRITELN <*Time infeasible*); 
tried := TRUE; 

IF excl_arc <> C3 
THEN BEGIN 


exci_arc := C3; 
goto ^7^; 

END; 

END; 

WRITELN < * s_nod * , s_nod : > , ’ d_nod * , d_nod : 5 ) ; 
WRITE (’Short aray * > ; 

FOR 11 : = no DOWNTO 1 DO 

WRITE (arcs_arayC i i 3 : 4) ; 

WRITELN; 

C We check the residual capacities of the arcs, 
and time feasibility before sending the fio > 

t,ravl_tim := 0; 

FOR 11 : = no DOWNTO 1 DO 

WITH arcCarcs_aray C 1 1 3 3 DO 
BEGIN 


^ 80 : 


tim_st := rk+travi_t im; 
passed := TRUE; 

FOR tt ;= round ( t im__st > TO round( t im_s t+d i s t > DO 
BEGIN 

IF res_capCtt3 <= 0 
THEN BEGIN 

1 1 m_s t : = 1 1 + 1 ; 
passed i == FALSE; 

END 

ELSE flo_sent := smai lerCf io_sent , res^capCttJ ) 

END; 

IF NOT passed 

THEN exci_arc ;= excl_arc + Carcs_aray C i i 3 3 ; 
travl_tim := (tim_st-rk> + dist; 

IF NOT passed AND ( ( trav i_t im<ava i l_t im> OR tried) 


THEN GOTO ^80; 
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AND NOT tried 


END; 


slack_tim := avail_tim - travl_tim; 

IF (slack__tim < 0> AND (mo 1 1 on_t i m<ava i l_t im> 

THEN BEGIN 

WRITELN (’Time infeasible (temp)’); 
GOTO 57^; 

END; 

{ else if motion_tim > tim then delays > 


IF avaii_tim < travl_tim 
THEN BEGIN 

WRITE (’Time i nf eas i b i 1 i ty ’); 
delay := travl_tim - avail_tim; 

END 

ELSE .delay : = 0 ; 

WRITELN (’delay * , de 1 ay : ^ : 2 > ; 

Update_de lay_i nfo(delay,f lo_sent , de lay_f req , wtd_de lay , ptr ) 
FOR i i : = no DOWNTO 1 DO 
WITH arcCarcs^aray C I 1 ] 3 DO 
BEGIN 

FOR tt := round( t im_st ) TO round( 1 1 m_s t-fd i st ) DO 
BEGIN 

res^capCtt] := res_capCtt3 - flo_sent; 
ioadCtt3 := ioadCtt] t flo_sent; 

END; 

WRITELN (’arc no ’ , arcs_aray C i i 3 : 2 , 

’ Res capa ’ , res__capC 1 1 3 : 5 : 0 , ’ load ’ , ioadC 1 1 3 : 5 : 0 ) ; 

WRITE (*tim s t : * , t im_s t : > : 0 ) ; 

1 1 m_d t : = 1 1 m_s t + d i s t ; 

WRITELN (tim_dt :5:0> ; 

period__cost := period_cost + f lo__sent*cost_^w; 

END; 

WRITELN (’flo sent ’ , f lo_sent : 5 : 0) ; 

ass i gnC 1 store . j store 3 : ~ ass i gnC i store , j store 3 -flo^sent; 
spnC i store3 . sup := spnC i store3 . sup - flo_sent; 
time := spnC istore3 . t im; 

WITH nodeCs_nod3 DO 
BEGIN 


requ_avai C t ime3 := requ_avai C t ime 3 - flo_sent; 
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cum_dem_sup_w cum__dem_sup_w - fio_sent; 

END; 

dmnC j store 3 . dem := dmn[ ) store 3 . dem - flo_sent; 
time r= dmnC j s tore 3 . t im'; 

WITH nodeCd_nod3 00 
BEGIN 

requ_avai Ct ime3 := requ_ava i C 1 1 me 3 fio_sent; 

cum_dem_sup_w := cum_dem_sup_w + fio_sent; 

IF cum_dem_sup__w > -0.001 

THEN dem__nod_set := dem__nod_set - Cd_nod3; 

END; 

fio_sent > 0.001 
THEN BEGIN 

WRITE <out2,f io_sent:6:2, * ’>; 

FOR 11 : = no DOWNTO 1 DO 

WITH arcCarcs_aray C 1 1 3 3 DO 

WRITE (out2,nodel:8,t im_s t : 8 : 2 ) ; 

IF no <> 0 

THEN WITH arcCarcs_arayCl 33 DO 

WRITELN(out2,node2:8, < 1 1 m___st +d i s t > :8:2) 
ELSE WRITELN <out2 , s_nod : 8 , spnC i store 3 . t im: 6 , 
d_nod : 8 , dmnC jstore3.tim:6> ; 

ENDj 

cn : = cn +1 ; 

END; C while (dem_nod_set <> C3> and (cn<200> do> 

IF dem_nod_set <> C3 

THEN WRITELN ('Exiting while loop !! emp wag dem unsatisfied'); 
WRITELN (out? , per i od_cost : 1 2 ; 2 > ; 
total_cost ;= total_cost + period_cost; 

600; IF period_no* period_lenth < T_hor 
THEN BEGIN 

penod^no := period_no + 1; 

GOTO ?0; 

END; 

REWRITE (outl . 'outdelay.dat ') ; 

WRITELN (outl, 'Delay Frequency % Delay wtd % Delay*); 
cumsum := 0; tot_delay := 0; 

FOR 11 := 0 TO ptr DO 


BEGIN 
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cumsum := cumsum +■ da lay_f reqC i i 3 ; 
tot_deiay := tot_delay + wtd_de lay C i i 3 ; 

END; 

FOR 11 := 0 TO ptr 00 

WRITELN < out 1 , i 1 : 4 , * ’ , de iay_f reqC i i 3 : 4 . delay 

cumsum*! 00 : 1 2 : 2 , 


wtd^de lay C 1 i 3 / tot_de lay*! 00 : ! 2 : 2 > ; 
WRITELN (out!); 

WRITELN (out!, ’Total Frequency 
WRITELN (out!, ’Total Wtd delay 
WRITELN ( out !,’ Average Delay 
CLOSE (out!); 

CLOSE <out2); 

REWRITE ( out > , ’ out 1 oad . dat * > ; 
WRITELN <out3^,’Arc no Load 
FOR 11 : = 1 TO 2*no__arc DO 

WITH arcCi i3 DO 
BEGIN 


’ , cumsum : 4 > ; 

’ , t ot_de i ay : 6 : 2 > ; 

’ , tot_de lay /cumsum 


Time ’ ) ; 


max load := 0; 

FOR tt ;= 1 TO T_hor DO 
BEGIN 

IF maxioad < loadCtt3 
THEN BEGIN 

max load := loadCtt3; 
t ime : = t t ; 

END; 

END; 

WRITELN ( out 5 , 1 i : 6 ,maxl oad :8:0,time:8>; 

END; 

CLOSE (out5>; 

REWRITE<out4, ’ da i 1 y load . d ’ ) ; 
per i od_no : = ! ; 

WHILE period^^no * pen od^^lenth <= T_hor DO 
BEGIN 

WRITELN (out4 , * Per 1 od_no : ’ , per i od_no : 2 > 

WRITELN (out4,’Arc No. Time’); 

WRITE <out4, • ’ > ; 

FOR tt ;== (period_no-! )*per iod_lenth + 1 
pariod^no * period_lenth 


f req[ i i 3 / 


: 6 : 2 ) ; 


to 

DO 


WRITE <out4,tt:4>; 
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WRITELN (out4>: 

FOR 11 : = 1 TO 2 * no_arc DO 

WITH arcC i i 3 DO 

BEGIN 

WRITE (out4,ii:4,' ’); 

FOR tt := ( per 1 od_no-l > *per 1 od_i enth + 1 to 

period_no * period_lenth DO 
WRITE (out4. loadCtt3:4:l) ; 

WRITELN (out4>; 

END; 

period_no := period_no + 1; 

END; 

CLOSE (out4> ; 

WRITELN (out^); 

WRITELN <out5, 'Total Cost : ' . tota l_cost : 1 2 ; 2 > ; 

CLOSE <out5); 

END. 



APPENDIX C 


PROBLEMS SOLVED AND THEIR RESULTS 

Here we present results of the 6 more problems. All the 
problems solved, took a CPU time in the range of 27-32 seconds of 
HP 9000. The results obtained are more less similar to those 
presented in the chapter 3. However here we have noticed that for 
the problems having tighter capacities of the arcs have not only 
more delays but also more costs of transportation, which appears 
clearly in graphs. 


PROBLEMS 3 AND 4 

The problems 1 & Z are solved as 9 (instead of 4> period problems. 
Time horizon = 96 time units. Period Length = 20 time units 
FINAL RESULTS FOR THE PROBLEMS 3 AND 4 - 

Table C. 1 Final Results for the Problem 3 and 4 


Problem 3 


Empty Uagons 

Transportation Costs 

Delays 

128 . 

36,446 

4,26 

133 

53,968 

3.97 

178 

33,363 

3.83 

203 

53,103 

3.66 

128 

39,677 

4.43 

133 

38.370 

4.12 

178 

58.261 

5.94 

203 

37,347 

3.73 


Problem 4 






0 


INPUT FOR THE PROBLEMS 5 AND 6 
Table C. 2. 1 Supply at various nodes 


Supply Node 

Supply, time (day) 

5 Kanpur 

4.5; 6.4; 7,5; 5.6; 10,7; 5.9; 4.10; 

5 Calcutta 

6.1; 7.4; 12.5; 5.6; 2.11; 8,12; 

7 Vishakhap. 

8.2; 5.5; 4.2; 6.5; 7,6; 5.9; 6.10; 

8 Madras 

5,1 ; 4.2; 8.7; 5.8; 4.11 ; 6.12; 

1 5 Ahmadabad 

2,1; 8,5; 5.7; 6.8; 6,11; 7,12; 

1 4 Manmad 

5,1; 4,2; 5,5; 4.6; 5.9; 4.10; 

16 Nagpur 

8.4; 5,1; 5,7; 4.8; 7.9; 5.11; 10,12 

19 Bhopal 

7.1; 5.2; 10.5; 4.5; 6.6; 8.9; 5.10; 

20 Indore 

5.2; 6,1; 2.7; 8.8; 8,9; 5,10; 4,11; | 


Table C. 2. 2 Demands at various nodes 


Demand node 

Demand, time (day) 

2 Jhansi 

4.5; 6.4; 4.5; 6.6; 8.9; 5,10; 

18 Hydrabad 

5,2; 6.1; 5.7; 6.8; 6,11; 7.12; 5.9; 6,10 

4 Lucknow 

8.1; 7,4; 8.5; 6.5; 7,6; 8.7; 5.8; 4,11; 


6,12; 5.9; 4,10; 

9 Bangalore 

8.4; 5,1; 4.2; 5.5; 4.6; 7.9; 5.11; 10,12 

6 Jamshedpur 

5.1; 4.2; 5.7; 4,8; 2,11; 8.12 

1 Delhi 

5,1; 4.2; 8.5; 5.6; 4,5; 2,7; 8.8; 8,9 


5.10; 4,11 


1 2 Bombay 


8.2; 1^.5; 7.?; 5.6; 10,7; 5,9; 4.10 














FINAL RESULTS FOR PROBLEMS 5 AND 6 

Table C*2.4 Final Results for the Problem 5 


Empty Wagons 

Transportation Costs 

Delays 

205 

22,471 

4.56 

250 

21 , 568 

00 

• 

255 

21 ,042 

4.15 

280 

21,042 

4.17 


Table C.2.5 Final Results for the Problem 6 


Empty Wagons 

Transportation Costs 

Delays 

205 

22,595 

5.64 

250 

21,480 

5.64 

255 

21,155 

5.55 

280 

21,155 

5.52 
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Transport. Cost Vs. Wagons Available 



Ka of Wtfom ATallabls 

FIG. C.3 

Delays Vs. No. of Wagons Available 
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FIG. C.4 
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INPUT FOR THE PROBLEMS 5 AND 6 
Table C. 2. 1 Supply at various nodes 


Supply Node 

Supply, time <day) 

3 Kanpur 

4.3; 6,4; 7.5; 5.6; 10.7; 5,9; 4.10; 

5 Calcutta 

6.1; 7,4; 12,5; 5.6; 2,11; 8,12; 

7 Vishakhap. 

8.2; 5.5; 4.2; 6.5; 7,6; 5.9; 6.10; 

8 Madras 

5.1; 4.2; 8,7; 5.8; 4.11; 6.12; 

1 5 Ahmadabad 

2,1; 8,5; 5,7; 6,8; 6,11 ; 7,12; 

1 4 Manmad 

5,1; 4,2; 3.5; 4.6; 5.9; 4.10; 

16 Nagpur 

8.4; 5.1; 5,7; 4.8; 7.9; 5.11; 10,12 

19 Bhopal 

7.1; 5.2; 10,5; 4,5; 6,6; 8.9; 5.10; 

20 Indore 

5.2; 6,1; 2,7; 8,8; 8.9; 5.10; 4,11; 


Table C* 2. 2 Demands at various nodes 


Demand node 

Demand, time (day) 

2 Jhans i 

4.3; 6.4; 4,5; 6.6; 8.9; 5.10; 

16 Hydrabad 

5,2; 6,1; 5.7; 6.8; 6,11; 7,12; 5,9; 6,10 

4 Lucknow 

8,1; 7,4; 8,3; 6.5; 7,6; 8,7; 5.8; 4,11; 


6,12; 3,9; 4,10; 

9 Bangalore 

8,4; 5,1; 4.2; 5.5; 4,6; 7.9; 5,11; 10,12 

6 Jamshedpur 

5,1; 4.2; 5,7; 4,8; 2,11 ; 8,12 

1 Delhi 

5.1; 4,2; 8,5; 5.6; 4.5; 2.7; 8,8; 8.9 


1 2 Bombay 


5.10; 4.11 

8.2; 15.5; 7.5; 5.6; 10.7; 5.9; 4.10 












FINAL RESULTS FOR PROBLEMS 5 AND 6 

Table C.2.4r Final Results for the Problem > 



Table C. 2*5 Final Results for the Problem 6 


Empty Nagons 

Transportation Costs 

Delays 

205 

22,395 

5.64 

250 

21,480 

5.64 

255 

21,155 

5.35 

2B0 

21,153 

5.32 
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INPUT FOR THE PROBLEMS 7 AND 8 

Table C. 3*1 Supplies at various nodes 


Supply Node 

Supply, time (day) 

1? Bhusaval 

6.l! 7.2? 1,6; 4.7; 6.11; 7,12 

20 Indore 

8.?; ?.4; 8.8; 7.9; 9.10; 8,12; ?.l? 

? Kanpur 

4.?; ?,1; 4.6; ?.7; 4.8; 1.14; 4.1? 

4 Lucknow 

4,2; 2.3; 4,9; 2,10; 8,11; 7.12; ?,13 

9 Bangalore 

4.4; ?.?; 4.1; 4.6; ?.7; 4.14; 3.16; 4,11 

6 Jamshedpur 

8.2; 7.3; ?,4; 8.8; 3.9; 4.12; 2,13 

1 Delhi 

1,3; 4.1; 6.10; 7.6? 4.14; 3.13 


Table C. 3*2 Demands at various node 


Demand Node 

Demand, time (day) 

19 Bhopal 

6.2; 4.10; 4,6; 3.13; 7,11 

1 3 Ahmdabad 

1.3; 4,1; 3,8; 3.9; 4,13; 4,14 

7 Vishakha. 

6,4; 1,10; 2.6; 3,7; 3.11; 3.12 

1 1 Pune 

1.2; 2.3; 3.7; 4,8; 1,9; 2,14; 6.13 

17 Balharsha 

7,3; 2.1; 3,10; 3,6; 8,1? 

1 0 Daund 

1.4; 8,9; 3.11; 3,12 

1 4 Manmad 

3.2; 3.3; 2,7; 6.8; 3,13; 1.13 

16 Nagpur 

4.3; 3.1; 3,10; 3.6; 1.13; 2.11; 3,12 

2 Jhans i 

1,3; 7.4; 4,8; 4,9; 3,13; 9,14 

1 8 Hydrabad 

3.1; 4,2; 3,6; 7,7; 4,11; 4,12 


The time horizon = 1? days and period length = ? days. 
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Table C. 3. 3 Table for capacities of the arcs. 


Arc 

No . 

B 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

Pr. 

7 

B 

5 

25 

8 

20 

24 

5 

21 

6 

18 

18 

7 

6 

23 

Pr. 

8 


5 

20 

8 

20 

20 

5 

20 

6 

15 

10 

6 

7 

22 



FINAL RESULTS FOR THE PROBLEMS 7 AND 8 
Table C.3.4 Final Results for problem 7 


Empty Wagons 

Transportation Costs 

Delay 

157 

21,753 

5.25 

182 

20,790 

5.03 

207 

19,786 

4.61 

232 

17.843 

4.62 


Table C.3. 4 Final Results for problem 8 


Empty Wagons 

Transportation Costs 

Delay 

157 

22,187 

6.06 

182 

21,241 

5.63 

207 

20.223 

5.32 

232 

18.300 

4.72 























DISCUSSION 


Here we give minutes of discussion during the thesis presentation. 

Q. What IS the smallest unit of time for which scheduling is 

done ? What is the basis of deciding the its level 7 

Ans. We have used 1 unit time = 4 hrs. in our scheduling problem. 

The criteria for deciding the smallest unit of time are 

1. Travelling time for arc with smallest dist 2. Smaller the unit 

time finer the level of scheduling and more computational time. 

Q. Where do we use simulation in the methodology 7 

Ans. During routing we try incorporate to have essence of real 

life logic and relationships. Further, we perform experiments we 

the model developed so as to decide what fraction of capacity of 

arcs per period should be taken as the capacity of arc per unit 

time. 

Q. What IS alternative methodology that can be used 7 
Ans. One can try to use Lagrangian Relaxation by relaxing some 
constraints, dualizing them and exploit some standard structure 
out of the problem. Also some other ways of dividing the problem 
into smaller subproblems, methods like Benders’ Decomposition may 
be tried. Also one can have analogy between our problem and job 
shop scheduling problem. Wagons are analogous to the jobs. 
Wagons moving over arcs is similar to the jobs loaded on the 



machines. Travelling time being analogous to processing time of 
jobs. Delivery dates for the goods are similar to due dates of 
the jobs. However it should be noted that scheduling many jobs on 
more machines (more than 3> is itself a tough problem to solve. 

Q. How the costs are expressed? Can they be expressed in terms 
of cost of moving a wagon per arc? 

Ans* Costs are expressed as cost of moving a train of wagons per 
arc. It can be expressed in terms of wagons also for going into 
finer details of scheduling. It will involve some component of 
cost which is fixed and some component which is proportional to 
the number of wagons in tram. 

ft. Is it necessary to take cost proportional to the distance 7 
Ans. It is not necessary to take cost proportional to the 
distance. In real world also cost may vary not only with distance 
but also depending on whether the route is electrified or uses 
diesel or coal engines. The algorithm will work such input also 
without any further modification. 



